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ABSTRACT 


We present new, more precise measurements of the mass and distance of our Galaxy's 
central supermassive black hole, Sgr A*. These results stem from a new analysis that more 
than doubles the time baseline for astrometry of faint stars orbiting Sgr A*, combining two 
decades of speckle imaging and adaptive optics data. Specifically, we improve our analysis 
of the speckle images by using information about a star's orbit from the deep adaptive optics 
data (2005 - 2013) to inform the search for the star in the speckle years (1995 - 2005). When 
this new analysis technique is combined with the first complete re-reduction of Keck Galactic 
Center speckle images using speckle holography, we are able to track the short-period star SO- 
38 (K-band magnitude = 17, orbital period = 19 years) through the speckle years. We use the 
kinematic measurements from speckle holography and adaptive optics to estimate the orbits of 
50-38 and 50-2 and thereby improve our constraints of the mass (Мл) and distance (Ro) of Sgr 
АХ: My, = 4.02 £0.16 £0.04 x 10° Мо and 7.86 +0.14 +0.04 kpc. The uncertainties in Mpn 
and R, as determined by the combined orbital fit of 50-2 апа 50-38 are improved by a factor of 
2 and 2.5, respectively, compared to an orbital fit of SO-2 alone and a factor of ~2.5 compared 
to previous results from stellar orbits. This analysis also limits the extended dark mass within 
0.01 pc to less than 0.13 х 109 Mo at 99.7% confidence, a factor of 3 lower compared to prior 


work. 


Subject headings: Galaxy: center — astrometry — techniques: high angular resolution — 
infrared:stars — galaxies: supermassive black holes — Galaxy: fundamental parameters 


1. Introduction 


Following the motions of stars in the center of our Galaxy has given many insights into the properties 
of the gravitational potential in which they move. The measurement of the high proper motions, and later, 
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accelerations of these stars implies that they move in the gravitational potential of а concentrated dark mass 
[Eckart et al2002). With further observations, 
these stellar motions have provided strong evidence for the presence of a supermassive black hole (SMBH) 
at the Galactic Center (Sgr АХ) with a mass of about 4 x 106 Mo. Once the star S0-2 went through closest 


approach in 2002, it was possible to fit its motion with a Keplerian orbit (Schódel et al.|2002; 
2003). In addition to the mass of the SMBH, stellar orbits with measured radial velocities have been used 


to determine the distance to the Galactic Center (Ro; 2003 2003). With SO- 


2's short orbital period, this star provides the best constraint on the mass of the central black hole and Ro 
from stellar orbits to date (e.g., (Gillessen et al.|2009b). The focus of recent work has 
been to continue assessing the central black hole's properties as well as exploring the potential for using the 
measurement of stellar motions to test general relativity. 


The mass of Sgr A* (Мь;) and the distance to the Galactic Center (Rọ) are both important ways of 
characterizing this unique region of our Galaxy and putting it in context with others galaxies. Measuring 
the mass of the central supermassive black hole allows it to be compared to supermassive black holes in the 
centers of other galaxies. With the mass of Sgr A*, the Milky Way can be added to observed correlations 
between the mass of the central SMBH and other galactic properties, such as velocity dispersion of stars and 
bulge luminosity see Kormendy & HoP0T3]for a review) 
Ro is a key parameter characterizing our galaxy's size, mass and kinematics. The adopted value of R, affects 
estimates of the Milky Way's rotation curve and thereby also measurements of the mass and shape of the 
dark matter distribution (e.g., [Olling & Merrifield|2000). An independent, very accurate measurement of 
Ко could possibly be used to calibrate stellar distance indicators, such as RR Lyrae and Cepheids, which 
are important steps оп the cosmic distance ladder (see[Reid|1993). А, additionally serves to calibrate the 
extinction towards the Galactic Center (e.g.,|Schödel et al.[2010). The mass-to-distance ratio of Sgr А" as 
derived from stellar orbits is also necessary to determine the predicted size of the black hole shadow that 
will be observed by the upcoming Event Horizon Telescope, which can be used as a null hypothesis test of 
general relativity (see [Psaltis et al.[2015). Finally, future tests of general relativistic effects on the motion 
of S0-2 depend on accurate measurements of the gravitational potential due to the SMBH in the Newtonian 
regime. 


Until now, our group has used only S0-2 to constrain R, and the mass of Sgr A*. This is because S0-2 
is unique: it is bright (14.2 in the K band) and has a short orbital period (16.2 years). We therefore have 
been able to track its motion since Keck observations of the Galactic Center began in 1995 such that our 
observations now cover more than one full orbit of this star. We would like to also use other short-period 
stars in the Galactic Center to help determine the gravitational potential, but ideally only those stars with 
high orbital phase coverage like 50-2. It has been shown for visual binary stars that if less than 40 - 5096 
of a body's orbit is covered by astrometric observations, the orbital parameters estimated from the data are 
systematically biased from their true values (Lucy|2014). We therefore only use additional short-period 
stars to constrain the black hole mass and К, that conservatively have at least this minimum orbital phase 
coverage. 


Achieving this minimum orbital coverage for other short-period stars is a challenge because those stars 
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are fainter than 50-2 by more than an order of magnitude, making them difficult to track in the first 10 years 
of speckle imaging data. In this work, we present a method to increase the orbital phase coverage of these 
faint stars at the Galactic Center through a complete reanalysis of this data set. In all past analyses of stellar 
orbital parameters using the speckle data, the data from each observation run has been treated independently. 
Stars are blindly searched for in the summed image from a given epoch of observation, as if the Galactic 
Center had never been observed before. No information from other observations is used in this search. In 
this work, we present a new methodology for analyzing the speckle images that does use information from 
the much deeper adaptive optics images and the vast improvement that has been made in the knowledge of 
the central black hole's properties. We use constraints on a star's orbit from the deep adaptive optics data to 
inform the search for the star in the earlier speckle years. 


As à pilot study for this new methodology, we apply this technique to 50-38, one of the three stars at 
the Galactic Center with an orbital period of less than 20 years (in addition to 50-2 and 50-102; see [Meyer] 
et al.[2012), with the ultimate goal of using this star as an additional constraint on the black hole mass and 
Ко. At a magnitude of К- 17.0, 50-38 is consistently detected in our deep adaptive optics data taken from 
2005 - 2013 but it has not previously been detected in our speckle imaging data taken from 1995 - 2005. 
50-38 is an ideal star for the application of this methodology because it is consistently detected in all 21 
adaptive optics images and its radial velocity has been measured (Gillessen et al.|2009b|and this work), so 
its orbit is well known even with just over 40% of its orbit covered by AO observations. Our results are also 
made possible by a new reduction of the speckle data using the more sophisticated reconstruction algorithm 


called speckle holography (Schódel et al.|2013). This is the first work that includes the speckle holography 
re-reduction of all Keck Galactic Center speckle data. 


The paper is organized as follows: Section D] describes the data sets used in this work, including the 
results of the new speckle holography reduction on the full set of Keck Galactic Center speckle imaging data. 
Section 3] describes the new methodology of analyzing the speckle images as applied to 50-38. Section [4] 
contains the results of the 50-38 orbital analysis, including improved constraints on mass and R, as well as 
extended mass. 


2. Data Sets 


This paper is based on three different types of data sets, which are briefly summarized here. The 
first type is previously reported speckle imaging data that, in this work, are re-reduced with the speckle 
holography technique (Schédel et al.|2013). The result is a higher-quality, deeper final image for each 
observation epoch (Section [2.1] and Table 2. The second is new and previously reported adaptive optics 
(AO) astrometry data (Section 2.2.1. New AO imaging observations (Table |2) are analyzed by methods 
used in prior work by our group to produce star lists of relative astrometric measurements. Initial star lists 
of relative astrometric measurements from re-analyzed speckle data and star lists from new and existing 
AO data are transformed to an absolute reference frame using updated measurements for the astrometric 
standards (see Appendix [A). The third data type is new and previously reported AO spectroscopic data. 


de 


These data are analyzed by standard techniques to extract radial velocity measurements (Section 2.2.2}, 
which are combined with radial velocity measurements from the literature. No astrometric measurements 
from the literature are used due to the difficulty in consistently transforming these measurements to our 
reference frame. These data and initial data reduction steps allow us to extend the time baseline for S0-38 
by a factor of two through analysis steps that are described in Section B] 


2.1. Speckle Imaging 
2.1.1. Existing Observations 


The speckle data sets used for this study provide astrometric measurements of the central ^ 5"x5" of 
the Milky Way between 1995 and 2005. While the details of these observations can be found in earlier 


papers from our group (Ghez et al.|1998 2000 2005aj [Lu et al.|2005] [Rafelski 
let аї. 2007), we provide a brief summary. Individual K-band (2.2 шт) frames were obtained with МІКС 


with a pixel scale of 20 milliarcseconds on the W. M. 
Keck I telescope. During each epoch of observation, between 2,000 and 20,000 frames were obtained using 
very short exposure times (0.1 sec) to freeze the distorting effects of the Earth's atmosphere. We begin 
our analysis of these data with individual frames that have had the instrumental effects removed (i.e., sky 
subtracted, flat fielded, bad pixel corrected, distortion corrected). Table [I] provides the dates and number of 
frames for each of the 27 epochs of observations. 


2.1.2. New Image Reconstruction and Initial Star Lists 


The individual short-exposure frames are combined via post-processing techniques that compensate for 
the blurring effects of the Earth's atmosphere. The result is a diffraction-limited image whose final Strehl 
ratio and depth depend on the algorithm used to make the combined image. Originally, our speckle images 


were reduced using the shift-and-add algorithm (SAA; 1998 2000 
2002 2005a} [Lu et al./2005} [Rafelski et al./2007). In this technique, the individual frames are 


shifted so that the brightest pixel of each frame is at the same position before the highest quality frames are 
averaged together (see |Christou[1991). With SAA, only the brightest speckle from each frame contributes 
to the diffraction-limited core of the final point spread function (PSF) and all other speckles are part of a 
substantial halo. In the present work, we employ a more sophisticated speckle reconstruction technique 


called speckle holography 1973} 1990 1998), as implemented by 
Schödel et al.|(2013). This technique involves the deconvolution of the observed, distorted images with the 
instantaneous PSF as measured from a set of reference sources. In this approach, the information from all 


the speckles contributes to the final diffraction-limited core. In the implementation of speckle holography 
presented here, the diffraction-limited core of the final PSF contains — 6.5 times more flux than the case of 
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the SAA PSFs (Strehl ratio in speckle holography ~ 0.4 versus Strehl ratio in SAA ~ 0.067] and a much 
greater fraction of the data obtained is used in the analysis (average fraction of frames used in speckle 
holography ~ 92% versus average fraction of frames used in SAA ~ 17%). See Figure[1]for a comparison 
between SAA, speckle holography, and AO images. For each epoch, both a final image and three subset 
images of similar quality are constructed. An early version of our analysis of these data was reported in 


our recent work on S0-102 (Meyer et al./2012). A subset of the epochs presented in that early analysis 
used a preliminary version of the algorithm described in Schódel et al.|(2013), with a smaller field of view 
and a different treatment of the PSF reference sources. In the present work, the full implementation of the 
speckle holography algorithm from [Schódel et al.| was used. Image quality was further improved 
by rebinning the speckle frames from the original 20 milliarcsecond pixel scale down to 10 milliarcseconds 
before applying speckle holography as was done in the original SAA analysis. This is the first complete 
application of speckle holography to the Galactic Center speckle imaging data set from Keck Observatory. 


After the final speckle holography images are constructed for each observation epoch, an initial search 
for stars in these images is performed. In this preliminary analysis, each observation epoch is treated inde- 
pendently. The positions and fluxes of stars in the field are determined by the PSF fitting program StarFinder 
(Diolaiti et al.[2000). In the running of StarFinder to detect sources in both the main and subset images, we 
choose to set the correlation threshold to 0.8 instead of the value of 0.7 used in Meyer et al.|(2012). This 
more conservative correlation threshold value is chosen to minimize the number of fake sources detected in 
this preliminary analysis of the images. This correlation threshold is also motivated by the fact that we will 
later go back and look for sources that this initial analysis missed (see Section В). To further ensure that 
each detection is robust and to estimate the measurement uncertainties, StarFinder is run in the same way 
on the three subset images. We keep a detection in our initial star list if it was detected in all three of these 
subset images as well the main image. The astrometry and photometry measurements are taken from the 
main image. The error on our astrometric and photometric measurements is taken as the error of the mean 
(standard deviation divided by v3) of the positions and fluxes determined by StarFinder in the three subset 
images. 


We use the initial star lists to determine the K-band limiting magnitude for "direct detections" in each 
speckle holography image. The K limiting magnitude is defined as the magnitude at which the cumulative 
distribution function of the observed K magnitudes of stars in the initial list reaches 95% of the total sample 
size. Figure P]shows the K-band limiting magnitude of the SAA and speckle holography images for each 
observation epoch, as well as the limiting magnitude of the AO images for comparison. In our data, the 
median SAA K limiting magnitude is 15.7 and the median speckle holography K limiting magnitude is 
16.4. For comparison, the median K limiting magnitude of our AO data is 19.4. Table [I] summarizes and 
compares the SAA and speckle holography analyses of the 27 epochs of speckle data used in this work. 


'We note that these Strehl ratios are computed from the post-processed speckle holography and shift-and-add images. These 
Strehl ratios differ from those of the raw frames, with the SAA ratios being a closer representation. 


SAA Adaptive Optics 


July 2004 July 2005 


Fig. 1.— Comparison of a speckle shift-and-add final image (left), a speckle holography final image (mid- 
dle), and an adaptive optics final image (right). The speckle holography and SAA images are made from the 
same observed frames from July 2004 and the AO image is from July 2005. Note that between these two 
observations dates the relative positions of the stars are significantly different. Each image is 1 arcsecond on 
a side. S0-38 is directly detected in both the AO and speckle holography images, but not in the shift-and-add 
image. Its detected position is shown with the black arrow. With a K-band magnitude of 17.0, 50-38 is easily 
detected in the deep AO images but it is at the direct detection limit for the speckle holography images. 


2.2. Adaptive Optics 
2.2.1. Imaging Data and Astrometric Reference Frame 


New high-resolution images of the central ~ 10" x 10" of our galaxy were taken in 2012 and 2013 
using the laser guide star adaptive optics system on the Keck 
II telescope. The images were taken with NIRC2 (PI: Keith Matthews) in the K'-band (2.1 шт) and have 
a plate scale of 9.952 milliarcseconds per pixel (Yelda et al.[2010). The AO observations of this field from 
2004 through 2011 have been previously reported 
Yelda et al.[2014). The new AO observations from the years 2012 and 2013 were 


collected in the same manner as earlier observations. In this setup, each frame was taken with an exposure 
time of 2.8 seconds with 10 coadds (280 times longer than the speckle imaging frames). Astrometry and 
photometry are extracted using the same techniques reported in our previous work. Uncertainties in these 
measurements are also determined by running the same analysis on 3 subset images constructed with 1/3 of 
the frames used in the deep image. Table ]summarizes the results from the new AO observations. Together 
with existing star lists from earlier Keck AO data, these 3 new star lists bring the total number of epochs 
of relative AO astrometry available for this study to 23 epochs, from which 21 AO detections of S0-38 are 
unambiguously made (see Figure ПУ. The average magnitude ої 50-38 as derived from these detections 
is 17.0 0.1. 50-38 is well above the direct detection limit of the AO data, but comparable to the direct 
detection limit of the speckle holography data (see Figure B). 
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Fig. 2.— The direct detection K-band limiting magnitude of the imaging data used in this work (from Table 
1). The limiting magnitude is defined as the magnitude at which the cumulative distribution function of 
the K magnitudes of the directly detected sources reaches 95% of the total sample size. The red squares 
show this limiting magnitude as derived for the original shift and add reduction of the speckle data and the 
blue circles show the improved depth that comes from the new speckle holography reduction of the same 
data. The purple triangles show the direct detection limiting magnitudes of the AO images for comparison. 
The dashed line shows the average magnitude of S0-38 derived from the AO data. This shows that SO- 
38 is consistently detected in the deep AO images but is near the limit of direct detection in the speckle 
holography images. Note that the limiting magnitudes of the first two AO epochs are similar to those of 
the speckle holography epochs. This is due to the small number of frames going into these AO images (10 
frames used versus the typical 100-200 frames). This fact also explains why S0-38 is detected in 21 out of 
23 AO epochs. 


In addition to the central field, we also take observations designed to measure the astrometry in the near 
infrared (NIR) of a set of 7 SiO masers. We tie the NIR measurements of these masers to astrometric and 
velocity measurements made in the radio to construct an absolute reference frame with Sgr A* at rest at the 
origin (see Reid et aL 2007, Yeda erar рота, 
Here we repeat the exact procedure detailed in[Yelda et al.|(2014) to construct the absolute reference frame, 
adding three new epochs ої МІК maser observations from 2011 to 2013. These new maser observations 
are summarized in Table [5|in Appendix [А] Appendix [A]also provides the resulting absolute astrometry for 
the IR secondary standards that are used to combine all the speckle holography and AO star lists of relative 
positions together into a common absolute reference frame. 
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2.2.2. Spectroscopic Data 


New spectroscopic observations for S0-38 were obtained using the integral field spectrograph OSIRIS 
with the laser guide star AO system on the Keck 1 telescope on 2013 May 11-13. 
The 50-38 data were taken with the K broadband filter (Kbb; 1.965 - 2.382 шт) in the 35 milliarcsecond 
(mas) plate scale. This observational setup is designed to measure the CO bandheads (2.3,um) of the short- 
period, late-type star 50-38. The resulting spectrum has a resolution of К ~ 3,600. Thirty four 900-second 
exposures were taken, each dithered by a small amount (^ 0.2 arcseconds). In 6 of these frames, the AO 
correction was such that 50-38 is confused with a neighboring source ~ 100 milliarcseconds away. The 
remaining 28 frames are used to obtain the spectrum of S0-38, providing a total on-source time of 7 hours. 
For calibration purposes, we observed skies with the same exposure time (900 s). The AO performance was 
excellent during these observations, providing a PSF with a full width at half maximum (FWHM) of ~ 70 
mas. 


The analysis used to extract spectra closely follows previous analyses (Do et аї (2013). For the new 
Kbb data, we extract S0-38's spectrum using a 35 milliarcsecond aperture. The resulting spectrum is cal- 
ibrated by subtracting an annulus around the extraction aperture with a width of 70 milliarcseconds. The 
calibrated and normalized spectrum is shown in Figureß] This spectrum of SO-38 has a signal-to-noise ratio 
of 10, computed over a small, featureless wavelength range (between 2.212 and 2.218 microns). The radial 
velocity (RV) of S0-38 is then determined by measuring the cross-correlation of the extracted spectrum and 
a template spectrum in the wavelength region around the CO bandhead absorption features (from 2.285 to 
2.340 jum). We use the spectrum of the M3II giant HD40239 observed by the SPEX telescope 
as the template spectrum. first reported the detection of the CO bandheads 
in this source, and the RV reported in that paper is also used for our orbital analysis. The observed radial 
velocity is finally transformed to the local standard of rest reference frame by adding a correction of 27 km 
sl. The error on the RV measurements is taken as the root mean square of the RV measurements made 
on three subset spectra that are each made with ~1/3 of the individual frames. From this we obtain an 
LSR-corrected RV of -111 + 25 km s^! for S0-38 in May 2013. 


Table 1. Summary of Speckle Imaging Observations 
Date Frames Frames Used Nstars Strehl Ratio Кіт (mag) Opos® (mas) 
(UT) (Decimal) Available Holo SAA Holo SAA Holo SAA Holo SAA Holo SAA Data Source“ 
1995 June 9-12 1995.439 7532 4201 1800 199 151 0.23 0.06 16.0 15.6 0.5 1.1 (1) 
1996 June 26-27 1996.485 5391 4287 865 136 71 0.23 0.03 15.7 14.4 1.9 1.8 (1) 
1997 May 14 1997.367 3400 3400 1837 229 139 0.47 0.04 16.2 15.5 0.8 1.3 (1) 
1998 April 2-3 1998.251 2744 2712 1639 170 83 0.31 0.04 16.0 14.8 1.3 1.5 (2) 
1998 Мау 14-15 1998.366 9708 9708 2102 231 126 041 0.04 16.6 15.6 0.5 1.3 (2) 
1998 July 3-5 1998.505 2352 2352 933 172 127 0.46 0.06 16.3 15.5 0.5 1.2 (2) 
1998 Aug 4-6 1998.590 19741 12061 1933 235 172 0.48 0.06 17.0 15.7 0.3 0.8 (2) 
1998 Oct 9 1998.771 2548 2166 1082 161 120 0.58 0.07 16.1 15.4 0.5 1.5 (2) 
1999 May 2-4 1999.333 9800 9775 1857 263 183 0.51 0.07 16.9 15.8 0.4 1.1 (2) 
1999 July 24-25 1999.559 5684 5684 2108 273 232 0.50 0.09 16.9 16.0 0.3 0.9 (2) 
2000 April 21 2000.305 2940 2940 805 122 63 0.18 0.04 15.5 14.3 1.4 2.1 (3) 
2000 May 19-20 2000.381 15680 15680 2492 293 242 0.47 0.08 17.2 16.0 0.3 0.7 (3) 
2000 July 19-20 2000.548 11172 10637 1581 269 194 0.34 0.07 16.6 15.7 0.7 1.1 (3) 
2000 Oct 18 2000.797 2352 2352 1517 140 TI 0.39 0.04 16.0 14.8 1.1 1.3 (3) 
2001 May 7-9 2001.351 7306 7306 1994 225 175 0.39 0.07 16.4 15.7 0.5 1.0 (3) 
2001 July 28-29 2001.572 6860 6565 1695 290 239 0.54 0.11 17.0 16.3 0.3 0.8 (3) 
2002 April 23-24 2002.309 13620 13501 1958 267 183 0.41 0.05 16.9 15.8 0.4 1.2 (3) 
2002 May 23-24 2002.391 18052 11800 1443 290 252 0.53 0.08 17.1 16.1 0.3 0.9 (3) 
2002 July 19-20 2002.547 8081 5518 1118 208 125 0.36 0.06 16.4 15.4 0.9 1.4 (3) 
2003 April 21-22 2003.303 9392 9392 1841 144 121 0.36 0.04 15.9 15.6 0.6 1.1 (3) 
2003 July 22-23 2003.554 5028 5028 1703 190 180 0.58 0.07 16.2 15.9 0.6 1.2 (3) 
2003 Sept 7-8 2003.682 6270 6228 1723 233 182 0.35 0.07 16.4 15.8 0.6 12 (3) 
2004 April 29-30 2004.327 9289 9289 1423 236 185 0.46 0.08 16.7 15.8 0.4 0.8 (4) 
2004 July 25-26 2004.564 13110 12920 2161 293 200 0.43 0.08 16.8 15.9 0.4 0.9 (4) 
2004 Aug 29 2004.660 2850 2850 1301 203 167 0.31 0.08 16.2 15.7 0.5 1.3 (4) 
2005 April 24-25 2005.312 10553 10553 1679 226 162 0.37 0.06 16.6 15.7 0.4 1.0 (5) 
2005 July 26-27 2005.566 6080 5683 1331 171 111 0.32 0.05 16.1 15.3 0.5 1.2 (5) 


* Kris is the magnitude at which the cumulative distribution function of the observed K magnitudes reaches 95% of the total sample size. 


bPositional error taken as error on the mean from the three sub-images in each epoch and includes stars with К < 15. 


*Data originally reported in (1)[Ghez et al.|[1998), (2) |Ghez et al.|[2000), (3)[Ghez et al.](20053), (4) [Lu et al.](2005}, and (5)[Rafelski et аї (2007 
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Fig. 3.— Left: Observed spectrum of 50-38 as well as the template spectrum (shifted to a radial velocity of 
0 km s^!; used for the cross correlation showing the wavelength range over which this 
analysis is performed. Dotted vertical lines indicate the rest wavelengths of the CO bandhead absorption 
features seen in both the template and in S0-38's spectrum. Right: The measured cross correlation coefficient 
between the template spectrum and S0-38's spectrum over a range of radial velocity shifts of S0-38. These 
velocity shifts are in the local standard of rest reference frame. 
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Table 2. Summary of New AO Imaging Observations 
Date Date Frames Frames FWHM Strehl Nstars Кіт? с РАИ 
(ОТ) (Decimal) Obtained Used (mas) (mag) (mas) 
2012 July 24 2012.562 224 162 58 0.29 2344 19.6 0.06 
2013 April 25-26 2013.318 268 140 67 0.17 1432 18.8 0.11 
2013 July 20 2013.550 239 193 58 0.29 2492 19.6 0.08 


* Kji, is the magnitude at which the cumulative distribution function of the observed K magnitudes reaches 95% of 
the total sample size. 


bPositional error taken as error on the mean from the three sub-images in each epoch and includes stars with К « 15 
and with an average distance « 4 arseconds from the black hole. 


—12-— 


3. Data Analysis 


The goal of the work presented here is to increase the time baseline of observations of S0-38 so that 
more than ~40% of its orbit is covered and its information can be combined with that ої 50-2 to get further 
constraints on the gravitational potential in the Galactic Center. While the procedure described in Section 
2 results in unambiguous detections of 50-38 in all epochs of deep AO data (2005-2013), the same is not 
true for all epochs of speckle holography data (1995-2005). We refer to the AO and speckle holography 
detections from Section 2 as "direct detections." In order to add more speckle holography detections of 
50-38, we use a new methodology that takes advantage of the knowledge of (1) the black hole properties 
from 50-2 and (2) S0-38's orbital motion from the AO epochs to search for 50-38 in the speckle holography 
maps. We refer to the detections that result from this procedure as "prior-assisted detections," in contrast to 
the "direct detections." 


The data analysis associated with establishing and using the prior-assisted detections of S0-38 is de- 
scribed in the following four sections. In Section [3.1] we discuss why S0-38 is chosen as the first star to 
which the new methodology is applied. In Section 3.2] we identify stars that are similar in magnitude to 
50-38 and that may have overlapped with 50-38 on the plane of the sky between 1995 and 2005. By tracking 
these stars along with S0-38, we can better determine to which star any new speckle holography detections 
should be assigned. In Section [3.3] we create positional predictions for 50-38, in addition to the stars in the 
sample, based on existing knowledge of the orbital parameters, and use these predictions to search for 50-38 
in the speckle epochs. In Section all the AO and speckle holography astrometry along with the RV 
measured for S0-38 are simultaneously fit with S0-2's measurements to further constrain the gravitational 
potential through which these stars move. 


3.1. Selecting First Star to Apply New Methodology 


In this work, we chose a single star as a pilot test for our new methodology of searching for new 
detections in speckle holography data. We limit the possible stars to those for which full orbital phase 
coverage may be acheived, including both turning points, by adding the 11-year speckle holography time 
baseline to the 9-year AO time baseline. This puts an upper limit of 20 years on the orbital period of 
candidate stars. The goal of this orbital period criteria is to select a star that would give the most information 
about the central gravitational potential once new speckle holography detections are included in its orbital 
fit. There are two stars not consistently directly detected in the speckle holography images that meet this 
constraint: 50-38 (P = 19 years) and 50-102 (P = 11.5 years; Meyer et al.[2012). We chose to apply the new 
methodology to 50-38 for the following reasons. First, 50-38 is consistently detected in every one of our 
deep adaptive optics data sets (23 astrometric measurements) while S0-102 is detected in just over half (14 
astrometric measurements). The reason that 50-102 is often lost in the adaptive optics data is that its orbit 
moves through areas that are very crowded by other stars and by the МІК component of Sgr А" itself. 50-38 
also moves through a similar area on the plane of the sky as it goes through closest approach, but it spends 
the majority of its orbit away from the very central area in which S0-102 moves, thanks to its longer orbital 
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period/semi-major axis and its projection on the plane of the sky. In the speckle years specifically, S0-38 
is predicted overall to be in a less crowded region of the sky than S0-102, which makes finding additional 
prior-assisted speckle holography detections more feasible. Additionally, the radial velocity of 50-38 has 
been successfully measured (this work and |Gillessen et al./2009b), whereas the radial velocity of S0-102 
has not. This is largely thanks to the CO band heads present in S0-38's late-type star spectrum. The distinct 
shape and depth of these features make them easier to detect than the Br-» line present in the spectra of 
early-type stars. In the specific case of the K broadband spectroscopic data used to measure S0-38's RV in 
this work (Section 2.2.3), we are пої able to measure the spectrum ої 50-102 because it is confused with 
nearby sources in most frames. S0-38 was near furthest approach and in a more isolated region at the time 
of its spectroscopic measurement, making the extraction of its spectrum possible. We therefore choose to 
apply the new methodology to 50-38. 


3.2. Identifying Sources Potentially Confused with S0-38 


We construct a sample of stars that may have overlapped with S0-38 on the plane of the sky between 
1995 and 2005. The orbits of these stars are tracked along with that of S0-38 to avoid mismatching a star 
with any new speckle holography detection. An initial sample set of possibly confusing sources is identified 
through a radius and magnitude cut, where the radius used here is with respect to Sgr A*. The radius is set 
by how far away a star could be in the AO epochs from 2005 - 2013 and still have traveled into the central 
0."3 (0.012 pc for Ro = 8 kpc), where S0-38's orbit lies, during the 1995-2005 speckle epochs. In order 
to set a radius cut, we assume a representative velocity for the other stars in this central region. We make 
a conservative estimate of this velocity by taking 3 times the measured velocity dispersion at the projected 
distance of 0."3 from the central black hole (c = 400 km/s; [Do et al.]2013). While this assumed velocity is 
not an exact calculation of the velocities of stars in this region, it is a conservative value for the velocities of 
possibly confusing stars. 


We then assume that the star is constantly moving at this velocity and that this velocity is fully in the 
plane of the sky in order to calculate the maximum radius of possible confusing stars. The maximum radius 
is then: 

Tmax = З x сКаї 0.73) x 10 yrs + 0."3 = 0."62 (1) 


This equation says that a star starting at a distance Fmax from Sgr А" and moving at a constant velocity of 
3 x с(0."3) towards the black hole on the plane of the sky would be at a distance of 0."3 from Sgr A* after 
10 years. Any star that is within Fmax at some point in the AO epochs could have been inside the 0."3-radius 
circle in which S0-38 is orbiting at some point in the speckle epochs. 


From the stars within this radius, we then eliminate those sources whose magnitude is not consistent 


with S0-38's within 3 sigma (Kmag = 17.00 + 0.24). Finally, we also remove sources whose velocities are 
not within 90? of pointing to Sgr A* on the plane of the sky. The final set of possible confusing sources is: 
50-104, 50-102, 50-103, 50-37, 50-40, 50-42, 50-45, 50-49, and 50-23. 
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3.3. Extending S0-38's Astrometric Orbital Phase Coverage using Orbital Priors 


We apply a new methodology of analyzing the speckle holography data to find sources not detected in 
the initial, blind analysis. Our technique uses the information in the adaptive optics astrometric and radial 
velocity measurements of S0-38 as a starting point of where to look for this star in the speckle holography 
images. The orbits of the sample of possibly confusing stars described in Section B.2]are tracked along with 
that of 50-38 so that any new speckle holography detection is not misidentified. The basic steps of this new 
methodology as applied to S0-38 are as follows: 


1. Perform an orbital fit ої S0-2's direct astrometric detections from AO and speckle holography epochs 
and radial velocities, with all 6 orbital parameters and all 7 gravitational potential parameters left free. 


2. Perform an orbital fit of the direct astrometric detections and radial velocity measurements (if avail- 
able) of 50-38 and the possibly confusing stars (Sectionß.3.1). In these fits, the gravitational potential 
is fixed to the best-fit solution from SO-2’s orbital fit. For the sample stars, both speckle holography 
and AO direct detections are included in this fit. For 50-38, only AO direct detections are included. 


3. From the probability distributions of the orbital parameters found in each orbital fit, determine the 
probability distribution of the (X, Y) (on-the-sky) position of each star in the next earlier speckle 
holography epoch relative to those epochs included in S0-38's orbital fit (i.e., start with the most 
recent speckle holography epoch and move backwards in time in subsequent steps). 


4. Search for а speckle holography detection of 50-38 in the image from this epoch. We first check 
whether S0-38 is detected in the initial analysis of the image. If it is not in the initial star list, then 
apply the newly-developed StarFinder Force 1 algorithm to the image (Sectionß.3.2). This algorithm 
takes as input the predicted pixel position of the star from the probability distribution of X and Y and 
searches the speckle holography image for a point source around that position that had not previously 
been detected in the initial star list. 


5. If anew detection is found that meets a minimum correlation (0.3) and signal-to-noise threshold (30) 
and it is not blended with another known star, add that point to the list of detections of 50-38 and refit 
the orbit with that point (Section B-3.1). If no detection is found (1.е., the PSF fit in Starfinder Force 
1 does not converge), then search again in the next earlier speckle holography epoch. Note that the 
potentially confusing stars are not searched for in the speckle holography images and their orbital fits 
are therefore not updated. The majority of these stars are typically directly detected in the speckle 
holography images because their orbital periods are longer than S0-38's and they move through less 
crowded regions on the sky. 


6. Repeat steps 3 - 5 until all speckle holography images are searched for S0-38. 


7. After all new prior-assisted detections are identified, the black hole parameter probability distributions 
are determined by a simultaneous fit with 50-2 (Section[.4). 
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3.3.1. 6-Dimensional Orbital Fits and Positional Priors 


Orbital fits are carried out several times in our analysis. The initial task in the new analysis is to fit the 
orbits of all stars in the sample. For each new speckle holography detection, the orbit of S0-38 is also fit 
again to improve the prediction for earlier speckle epochs. All the orbital fits performed in steps 2 - 5 assume 
that the star moves in the gravitational potential of a single point mass that is the central supermassive black 
hole (BH), Sgr A*. In all orbital fits used to inform the search for the star in speckle holography images, we 
use the information from 50-2 to determine the probability distributions of the 7 parameters that describe 
the gravitational potential: the black hole's mass (М), position on the sky (xo, yo), 3D velocity (V, Vy, 
VJ), and the line of sight distance (Ro). The details of how 50-2 is used to determine the BH parameter 
probability distributions as well as our new S0-2 measurements are presented in Appendix [B] We then use 
the measurements of the star's position and radial velocity along with the constraints from 50-2 on the black 
hole parameters to determine probability distributions of the star's 6 Keplerian orbital parameters: orbital 
period (Р), time of closest approach (7,), eccentricity (e), inclination (i), argument of periapse (w), and angle 
to the ascending node (0). 


Our approach of using externally constrained black hole parameters speeds up the process of determin- 
ing positional priors for 50-38 and the other sample stars by a factor of ~ 5, allowing these orbits to be fit 
in a few hours instead of 1 day. The external constraints are implemented by stepping through a grid in the 
7D BH parameter space with 5 steps in each dimension, placed at -1, -0.5, 0, 0.5, and 1 sigma from the 
mean in that dimensionP] Each of the 57 - 78,125 points in this 7D grid is a single combination of values of 
the 7 BH parameters. These combinations are assigned weights based on the probability of those values as 
determined by 50-28 The weights of these combinations are very unevenly distributed because some sets 
of BH parameters are highly correlated. In order to not include combinations of BH parameter values that 
have a very low probability due to these correlations, we keep only the combinations that together make up 
95% of the total cumulative weight. This weight cut reduces the number of individual combinations by a 
factor of ~14 to 5,500. 


We then perform a series of individual fits to the star's detections in astrometry and radial velocity, one 
for each BH parameter value combination that made it past the weight cut. In each of these individual fits to 
the star's data, the 7 BH parameters are fixed to the values in each combination and the 6 orbital parameters 
are left free. Each individual fit then results in a set of samples of the posterior distribution of the star's 
6 orbital parameters. These individual fits are performed using the Bayesian multimodal nested sampling 
algorithm MULTINEST [Feroz et al.[2009), which was implemented in our code by 
Leo Meyer. This algorithm samples the posterior distribution more efficiently than traditional Markov chain 
Monte Carlo (MCMC) sampling schemes, especially when the distribution is multimodal. The MULTINEST 
algorithm has the advantage over our previous orbit-fitting method (see in that the best- 


?We choose to use a 1 sigma range in this grid spacing to make the fitting process computationally feasible. 


?The 7D joint probability distribution of the BH parameters as determined by 50-2 is approximated by a multivariate gaussian. 
This approximation represents the distribution well (see Figurd8) 
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fit orbital parameters and their probability distributions are found simultaneously and with much less total 
computing time. When the MULTINEST is run, uniform priors are applied to every free parameter, with 
ranges chosen such that the full range of the posterior distribution is probed. The results of the individual 
fits are then combined in order to get the final probability distributions for the orbital parameters of the star. 
This is achieved by combining the sets of samples of the posterior distributions found in each individual 
fit with the weight assigned to that fit’s unique BH parameter value combination. Finally, the probability 
distributions of the individual orbital parameters are determined by constructing histograms of the weighted 
posterior samples. 


The result is a probability distribution for each of the 6 orbital parameters determined by the star's 
detections, with the BH parameter information determined only by S0-2's detections. From the final prob- 
ability distributions of all the BH and orbital parameters we can then derive the probability distribution of 
the (X, Y) position on the sky of the star in the speckle holography observation epochs. These probability 
distributions of X and Y constitute our prior knowledge about the position of this star in all earlier epochs 
and are used as inputs into the StarFinder Force 1 algorithm (Section [3.3.2]. 


Although orbital fits are performed for all the stars in the sample, there are only two that came close 
enough to be potentially confused with 50-38: 50-104 (in 2002) and 50-49 (in 1998). Figure M| shows the 
probability distributions of X and Y for S0-38 and these confusing sources in one of the epochs in which 
these stars are close on the plane of the sky (July 2002), as well as a more typical epoch where S0-38 is 
isolated from the possibly confusing sources (September 2003). Any new detection that is found in the 
images from these observation epochs is only used in the fit of 50-38 if both 50-38 and the nearby source 
are separately detected. 


The initial orbital fit of 50-38 includes only the AO direct astrometric detections and radial velocity 
measurements, while the orbital fits of the possibly confusing stars include both AO and speckle holography 
direct detections. S0-38 moved so quickly in the speckle years compared to other stars in the sample that its 
direct speckle holography detections are generally too far away to associate with 50-38 by eye. An initial fit 
of S0-38's orbit with the gravitational potential fixed to S0-2's best-fit solution gives a reduced chi-squared 
of 12.2, implying that the AO errors determined by the 3 subset image method are underestimated. In order 
to correct this, we add in quadrature to the 3 subset image errors an additive error term constant for all 
AO epochs so that the reduced chi-squared of the AO-only fit equaled 1.0. This additive term is set to 1.7 
mas, compared to an average 3 subset image error of 2.9 mas for S0-38's AO direct detections. This same 
procedure is not required for S0-2 because the reduced chi-squared of the orbital fit of its AO and speckle 
holography direct detections is already ~ 1. 


3.3.2. StarFinder Force I Algorithm 


We can now use the probability distributions of S0-38's position on the sky in a given speckle epoch to 
inform our search for the star in that image, while also tracking the orbits of the potentially confusing stars 
included in our sample. We first check whether there is a direct detection in the initial star list for that epoch 
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Fig. 4.— Probability distributions of X and Y for S0-38 and nearby sample stars in two example epochs: 
September 2003 (/eft) and July 2002 (right). The pixel position is shown on the image, which has a plate 
scale of 10 milliarcseconds per pixel. In September 2003, the prediction of the (X, Y) position of S0-38 is 
well-isolated from predicted positions of all the other stars in the sample. This level of isolation is typical 
throughout the speckle epochs. In July 2002, however, S0-38 is predicted to be within 70 milliarcseconds 
of S0-104. These two sources are not separately detected in the 2002 speckle epochs due to their close 
proximity, so no new detection is assigned to 50-38 in these epochs. 


that coincides within uncertainties with the predicted position of S0-38. If a detection is not found in the 
initial analysis, then the newly-developed StarFinder Force 1 algorithm is applied. 


The StarFinder Force 1 algorithm is designed to search for any additional point source that is not in the 
initial star list. The main inputs to the algorithm are the best guess whole pixel position of the source (Xo, 
Yo), the size of the search box over which to initially look for the source, and the positions and fluxes of the 
stars in the initial list for that epoch. (Xo, Yo) and the search box size are chosen so that the initial search is 
done over the whole 3 sigma range of the probability distribution of X and Y, but in the end new detections 
are allowed to be found outside the search box. 


StarFinder Force 1 uses many of the sub-routines of the original StarFinder code and follows the same 
basic steps, outlined here: 


1. First, the code takes a subset of the main image centered on the input position (Xo, Yo). The size of 
this subset depends on the full width at half maximum (FWHM) of the point spread function (PSF). 
For a typical PSF FWHM, the size of the subset image in this correlation step is 5 pixels on a side. 
It then finds the position at which the correlation between this small subset image and the inputted 
PSF is maximized. The search box input parameter sets the offsets over which the PSF and the subset 
image are cross correlated. 


2. Next, the code takes a larger subset of the image (for a typical PSF FWHM, this subset image is 
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9 pixels оп a side) and performs the initial fit of that subset image. The free parameters of this fit 
are the positions and fluxes of all stars within the subset image (including any sources found in the 
initial star list plus the additional source). The initial guess for the position of the additional source is 
the position of maximum correlation found in the previous step. The position of the additional point 
source is allowed to move away from this initial guess during the fit and is not constrained by the 
inputted search box or the previously predicted (X, Y) position. Any contributions of flux in the subset 
image from stars with positions outside the subset image are considered fixed and are subtracted from 
the subset image before the fit is performed. 


3. Finally, a series of 2 iterative fits is performed over a slightly larger subset image (for a typical PSF 
FWHM, this is 11 pixels on a side). If these iterative fits converge, then the best-fit position and flux 
of the newly-detected source are returned. To assign this new detection to 50-38, we require that the 
new detection (1) has a correlation above 0.3, (2) has a flux more than 3c above the noise, and that 
(3) S0-38 is not confused with any other source in the image. 


We generally cannot assign errors to the points detected with StarFinder Force 1 using the three subset 
image method that is used in the construction of the initial star lists. Since the magnitude of S0-38 is close to 
the detection limit in the deep speckle holography images of each epoch, it is not consistently detected in the 
three subset images that use only 1/3 of the frames. We therefore assign the error of each StarFinder Force 1 
astrometric point so that the reduced chi-squared of the orbital fit with the additional point is approximately 
1. This is done on an epoch-by-epoch basis, so the errors of each new point are set independently. In this 
way, an updated orbital fit with the new detection can be done before the search for the star in the next earlier 
speckle holography epoch is performed. If 50-38 is detected in the initial analysis of that holography epoch 
(and therefore is already detected in all three subset images), then the errors determined by the three subset 
image method are used in the subsequent orbital fit. 


3.4. Orbital Fits That Use S0-38 to Constrain the Gravitational Potential 


The result of applying the method outlined above to all speckle holography images is a final set of 
astrometric measurements of 50-38, including both the direct detections and the prior-assisted detections. 
With a higher orbital phase coverage that covers both closest and furthest approach, we now use S0-38 to 
constrain the gravitational potential. 


Before performing these orbital fits, we reassign the errors of S0-38's direct and prior-assisted speckle 
holography detections to a single value for all epochs. The single value is chosen so that the reduced 
chi-squared of an orbital fit of 50-38 with the black hole potential fixed to the best-fit solution from 50- 
2 equals 1. This reassignment is done so that the errors on the position of S0-38 are not underestimated 
in the epochs of the direct detections, whose errors were originally assigned by the three subset image 
method. The errors determined by the error of the mean of the measurements on three subset images are 
expected to be on-average slightly underestimated and are also a noisy estimator of the true error of the main 
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map measurement. The reassignment therefore ensures that the errors on S0-38's positions, as well as the 
resulting errors on the black hole parameters, are not underestimated. This issue of underestimated error 
bars is already addressed in the AO epochs by the inclusion of an additive error term so that the AO-only 
orbital fit gives a reduced chi-squared of 1. Additionally, with this reassignment the errors of the direct and 
prior-assisted detections are treated consistently in that both are set by requiring the reduced chi-squared to 
equal 1.0. 


Once S0-38's speckle holography errors have been reassigned, we perform two Keplerian orbital fits 
in which S0-38 is used to constrain the gravitational potential parameters: (1) a 13-dimensional fit of only 
the orbit ої S0-38 and (2) a 19-dimensional fit of the orbits of both 50-2 and 50-38. The 13 free parameters 
in the first fit are the 7 black hole potential parameters and the 6 Keplerian orbital parameters that describe 
the motion of S0-38. This fit is used to compare the gravitational potential solutions determined by 50-38 
alone with those determined by S0-2 alone. The 19 free parameters in the second fit are the 7 black hole 
potential parameters and a set of 6 Keplerian orbital parameters for each star. This simultaneous fit of the 
orbits of S0-2 and S0-38 is used to determine the final best-fit values and errors of the black hole potential 
parameters. Both the 13- and 19-dimensional fits are performed using MULTINEST. 


We finally perform a simultaneous fit of the orbits of 50-2 and 50-38 that includes deviations from a 
pure Keplerian orbit due to a distribution of extended dark mass. From this fit, we find an upper limit for the 
amount of extended mass within the orbits of these stars. Following [Ghez et al.|(2008), we set the extended 
mass distribution to follow a power-law density profile such that the total enclosed mass is given by: 


3-У 
M(« r) 2 Ми Мом (є ro) (2) (2) 


The power law slope, y, is fixed to 1.5 in the orbital fits presented here, though (2008) showed 
that the extended mass upper limit within the small radius enclosed by stellar orbits did not depend strongly 


on the value chosen for ~y in the range from 0.5 to 3. We set the outer radius cutoff of the extended mass 
distribution to 3.4 x 10!! km = 0.011 pc, such that it encloses the apoapse distances of 50-2 and 50-38 
within uncertainties. Our data is not sensitive to the mass outside of the orbits of these stars. For this fit, we 
also set the characteristic radius in the equation above, го», to the same value of 0.011 pc. We then use the 
extended mass parameters to find the extended dark mass within the apoapse ої SO-2’s orbit (0.01 pc), to 
compare with results from previous works. This orbital fit is performed using MULTINEST and has 20 free 
parameters: 7 black hole potential parameters, a set of 6 Keplerian orbital parameters for each star, and the 
total extended dark mass within г, (Мох (« r;)). We also perform an extended mass fit of the orbit of 50-2 
alone for comparison purposes, which is discussed in Appendix [B] 


As the final part of these analyses, we explore the additional uncertainties and biases arising from the 
construction of the absolute reference frame. In Appendix|C] we describe the details of a jack knife analysis 
performed on the 7 SiO masers used to tie our NIR observations to the absolute reference frame. The results 
of this jack knife analysis are used to correct the statistical bias (which is small compared to the orbital 
fitting uncertainties for Mp, and Ro) and add an additional source of error for all the black hole potential 
parameters. These bias shift and the additional errors are only significant for the parameters describing the 
position and velocity of Sgr A* on the plane of the sky. 
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3.5. Simulations to Investigate Statistical and Systematic Errors 


To confirm that the errors on the black hole parameters resulting from the simultaneous fit of S0-2 and 
50-38 are reasonable given our measurement uncertainties, we simulate mock 50-2 and 50-38 data. We 
then fit these mock data and compare the widths of the resulting Mj; and А, probability distributions to the 
results from the real data. 


For every speckle holography, AO, and RV observation epoch, a simulated measurement is created by 
drawing from a gaussian distribution. This gaussian distribution has an average equal to the best-fit model 
from the combined S0-2 and S0-38 orbital fit and a standard deviation equal to that epoch's measurement 
error. The mock data are then fit using MULTINEST, in the same way as the real 50-2 and S0-38 data. The 
widths of the resulting Мь; and R, probability distributions are quantified as the standard deviation of the 
posterior samples generated by the MULTINEST fit. One hundred mock data sets are created and fit in order 
to sample the possible widths of the Mpn and А, distributions well while not requiring large amounts of 
computation time. 


4. Results 
4.1. New Speckle Holography Detections of S0-38 


The application of the new methodology of analyzing speckle holography images on S0-38 results in 
12 new astrometric detections of this star. These new detections go back to May 1997, spanning nearly 
the entirety of our time baseline for speckle holography. Of these 12 new speckle holography detections, 
6 are direct detections found in the initial analysis of the speckle holography data and 6 are prior-assisted 
detections found using the force method. The 6 direct detections are in the 6 deepest epochs in which 
S0-38 was not confused with another source. All the prior-assisted detections have K magnitudes that 


are consistent with S0-38's average AO magnitude within +0.5, are 3 - 20 sigma above the noise, and 
have correlations ranging from 0.7 - 0.85 (well above the 0.3 minimum value). Adding these new speckle 
holography detections to the 21 adaptive optics detections of S0-38 brings the total number of astrometric 
measurements to 33 and almost doubles the time baseline that these measurement cover. Table[3]lists all the 
astrometry and radial velocity measurements of 50-38 used in this study. Appendix[B|lists the measurements 
of S0-2 used in the orbital modelling. With the new speckle holography detections, S0-38 has now been 
observed for over 80% of its 19-year orbit (see Figure [5). 


50-38 is not detected in the other 15 speckle epochs for a variety of reasons. 50-38 is confused with 
another star in 7 speckle epochs, either with a star of similar brightness (with 50-104 in April and May 2002 
and with 50-49 in July 1998) or with а much brighter star (with 50-1, K - 14.7, in October 2000 and with 
50-20, K - 15.8, in July 2005 and May and July 1999). 50-38 was not detected in 6 other speckle epochs 
because the PSF fit performed by Starfinder Force 1 did not converge (April 2005, September 2003, July 
2002, May 2001, April 2000, and April 1998). Finally, only 2 detections were rejected because they did 
not meet either the correlation or signal-to-noise threshold: June 1996 (signal-to-noise = 1.7) and June 1995 
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(correlation = 0.04). 
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Table 3. 50-38 Astrometric and Radial Velocity Measurements 
Date К! AR.A. ADec. AR.A. Error | ADec. Error у, V. Error | Reference? Detection 

(Decimal) (mag) (arcsec) (arcsec) (arcsec) (arcsec) (km s7!) (km s7!) Method 

1997.367 16.9 -0.229 -0.076 0.012 0.012 - - (1) Speckle (prior) 
1998.366 16.8 -0.188 -0.088 0.012 0.012 - - (1) Speckle (prior) 
1998.590 17.9 -0.188 -0.088 0.012 0.012 - - (1) Speckle (direct) 
1998.771 17.0 -0.164 -0.079 0.012 0.012 - - (1) Speckle (prior) 
2000.381 17.3 -0.125 -0.100 0.012 0.012 - - (1) Speckle (direct) 
2000.548 16.5 -0.113 -0.102 0.012 0.012 - - (1) Speckle (direct) 
2001.572 17.0 -0.068 -0.083 0.012 0.012 - - (1) Speckle (direct) 
2003.303 17.0 -0.007 0.018 0.012 0.012 - - (1) Speckle (prior) 
2003.554 16.8 -0.002 0.029 0.012 0.012 - - (1) Speckle (prior) 
2004.327 17.0 -0.067 0.060 0.012 0.012 - - (1) Speckle (direct) 
2004.564 17.0 -0.074 0.064 0.012 0.012 - - (1) Speckle (direct) 
2004.660 16.6 -0.075 0.086 0.012 0.012 - - (1) Speckle (prior) 
2005.580 17.1 -0.1285 0.0622 0.0045 0.0023 - - (1) AO 
2006.336 16.9 -0.1596 0.0616 0.0017 0.0017 - - (1) AO 
2006.470 17.0 -0.1647 0.0599 0.0020 0.0025 - - (1) AO 
2006.541 16.9 -0.1668 0.0593 0.0017 0.0017 - - (1) AO 
2007.374 17.2 -0.1933 0.0467 0.0017 0.0018 - - (1) AO 
2007.612 16.9 -0.1990 0.0437 0.0018 0.0017 - - (1) AO 
2008.260 - - - - - -185 70 (2) AO 
2008.371 16.9 -0.2140 0.0322 0.0019 0.0018 - - (1) AO 
2008.562 16.9 -0.2161 0.0311 0.0022 0.0017 - - (1) AO 
2009.340 17.1 -0.2361 0.0225 0.0020 0.0019 - - (1) AO 
2009.561 17.1 -0.2389 0.0185 0.0018 0.0018 - - (1) AO 
2009.689 17.1 -0.2426 0.0156 0.0019 0.0021 - - (1) AO 
2010.342 17.0 -0.2473 0.0090 0.0017 0.0017 - - (1) AO 
2010.511 16.9 -0.2461 0.0056 0.0024 0.0022 - - (1) AO 
2010.620 17.1 -0.2495 0.0045 0.0019 0.0019 - - (1) AO 
2011.401 17.0 -0.2585 -0.0119 0.0022 0.0018 - - (1) AO 
2011.543 16.9 -0.2576 -0.0127 0.0017 0.0017 - - (1) AO 
2011.642 17.0 -0.2558 -0.0136 0.0018 0.0018 - - (1) AO 
2012.371 17.0 -0.2545 -0.0260 0.0028 0.0020 - - (1) AO 
2012.562 17.0 -0.2570 -0.0285 0.0018 0.0017 - - (1) AO 
2013.318 16.7 -0.2523 -0.0391 0.0033 0.0020 - - (1) AO 
2013.360 - - - - - -111 25 (1) AO 
2013.550 16.8 -0.2555 -0.0413 0.0021 0.0020 - - (1) AO 


“Measurement originally reported in (1) this work and (2)|Gillessen et al.|(2009b). 


^V. values listed here are in the local standard of rest reference frame. 
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Fig. 5.— The best-fit orbit for 50-2 (blue line) and for 50-38 (red line) on the plane of the sky. These 
model orbit lines show the positions of these stars from 1995 to 2014. Both stars orbit clockwise on the 
plane of the sky. Closed circles indicate astrometric detections that were used in the orbital fits. Open 
circles indicate points that were not used in the fits because these astrometric detections are biased due to 
the proximity of other known sources on the plane of the sky. For S0-38, this consists of the two epochs 
of May and June 2002 in which the position of S0-38 shows a bias because it is blended with S0-104; in 
all other confusion epochs Starfinder Force 1 can only recover a position close to the neighboring source 
and 50-38 is undetected. Note in particular the paucity of detections around the point of closest approach 
at the bottom of SO-2’s orbit. The implication of the absence of useable detections in this part of S0-2's 
orbit is that 50-2 constrains the horizontal position of Sgr А" in our reference frame (хо) better than its 
vertical position (yo). S0-38's orbit on the plane of the sky is perpendicular to that ої S0-2, so it can provide 
additional constraints on yo. 


4.2. Gravitational Potential Parameters 
4.2.1. Values Derived from the Orbital Fit of 50-38 Alone 


We fit the orbit of 50-38 alone, leaving the orbital parameters as well as the black hole parameters free, 
to compare with the fit of SO-2’s orbit alone (discussed in Appendix [B). The motion of 50-38 alone is well 
fit by a Keplerian orbital model. Comparisons of the best-fit model from the fit SO-38’s data alone are shown 
in Figure |6|(middle and bottom rows). Similarly, the best fit to S0-2's data alone is displayed in Figure[7] 


-24- 


The orbital model fit of S0-38's motion alone provides independent estimates the black hole potential 
parameters Хо, Yo, Ух, and V, (see Figure 8) and Table 4). The other black hole parameters (V., Ro, Mon) 
rely on radial velocity measurement and are not well determined. This is because 50-38 has only two radial 
velocity measurements and these RV values are close to 0 km 8 !. These three parameters are still left free 
in this orbital fit of SO-38, but the uniform priors applied to these parameters are set to ranges that fully 
encompass at least the +40 uncertainties derived from the fit of 50-2 alone (see Table В. Those that are 


independently well fit are those that only require astrometric measurements (Хо, yo, Vx, Vy; see Figure[)). The 
value for these parameters preferred by 50-38 agree with the values given by 50-2 within 1 sigma, except for 
the V, parameter which differs by 2 sigma. Note that reference frame uncetainties are a significant additional 
error source for xo, Yo, Vy, and V, (see Figure [0), but are negligible for Мр and А, (see Appendix |C). 
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Table 4. Best-Fit Black Hole and Orbital Parameters as Derived From the Fit of S0-2 alone, S0-38 alone, 
and the Simultaneous Fit ої 50-2 and 50-38 


Best-Fit Parameter Values from Orbital Fits? 
Model Parameter (units) S0-2 only S0-38 only S0-2 and S0-38 


Black Hole Properties: 


Distance (kpc) 8.02 + 0.36 + 0.04 [6.5,9.5]P 7.86 +0.14 + 0.04 
Mass (1080) 4.12 +0.31 + 0.04 12.5, 5.51 4.02 £0.16 + 0.04 
X Position of Sgr A* (mas) 2.52+0.56+ 1.90 5.25 + 9.41 + 1.90 2.74 +0.50 1.90 
Y Position of Sgr A* (mas) 4.37 + 1.34 x: 1.23 6.85 + 5.00 + 1.23 5.06 + 0.60 + 1.23 
X Velocity (mas/yr) 0.02 + 0.03 + 0.13 0.40 + 0.70 + 0.13 0.04 + 0.03 + 0.13 
Y Velocity (mas/yr) 0.55 + 0.07 + 0.22 0.48 + 0.43 £0.22 0.51 +0.06 +0.22 
Z Velocity (km/sec) 154104 1-80, 4015 15.48 + 8.36 44.28 


S0-2 Properties: 


Period (yr) 15.90 + 0.04 15.92 + 0.04 
Time of Closest Approach (yr) 2002.343 + 0.008 2002.347 + 0.003 
Eccentricity 0.890 + 0.005 0.892 + 0.002 
Inclination (deg) 134.7+0.9 134.2+ 0.3 
Argument of Periapse (deg) 66.5+0.9 66.8 + 0.5 
Angle to the ascending node (deg) 227.9+0.8 228.0+ 0.5 


50-38 Properties: 


Регіоа (уг) 19.1+0.6 19.2 + 0.2 
Time of Closest Approach (yr) 2003.1 + 0.1 2003.19 + 0.04 
Eccentricity 0.80 + 0.02 0.810 + 0.004 
Inclination (deg) 170-4 1703 
Argument of Periapse (deg) 2030 12x21 
Angle to the ascending node (deg) 110+ 30 95 + 20 
Reduced Chi-Squared 1:5 0.9 1.3 


? AII best-fit values are shifted by the bias determined by the jack knife analysis of the reference frame (see Ap- 
pendix [C]. The first error term for each best fit value corresponds to the error determined by the orbital fit. For the 
black hole parameters, the second error term corresponds to jack knife uncertainty from the reference frame. 


‘These parameters are not well-determined by 50-38 alone, so instead of a best fit the uniform prior ranges that 
were used in the fit are reported. 
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Fig. 6.— Orbital fits of all speckle holography and AO astrometry and radial velocity measurements of 
the short-period star S0-38. The top row shows the measurements of S0-38's position and radial velocity 
along with the best-fit models of S0-38's orbit as determined by the fit of 50-38 alone (dash-dotted line) 
and the simultaneous fit of 50-38 and 50-2 (solid line). The middle and bottom rows show the residuals of 
the position and radial velocity measurements from the best-fit models. The middle row shows the residuals 
from the fit of 50-38 alone and the bottom rows shows the residuals from the simultaneous fit of 50-2 and 
50-38. The dashed lines in these rows indicate the 1-sigma uncertainties in the best fit models. In all plots, 
the detections of 50-38 are shown in different colors and shapes depending on the method in which they were 
detected: AO detections are purple circles, direct detections from the speckle holography images are teal 
squares, and prior-assisted detections from the speckle holography images are red triangles. The astrometry 
of 50-38 is plotted as the offset from the position of Sgr A*-radio, which is defined as the origin of our 
absolute reference frame. 
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Fig. 7.— Orbital fits of all speckle holography and AO astrometry and radial velocity measurements of the 
short-period star 50-2. The top row shows the measurements of S0-2's position and radial velocity along 
with the best-fit model of S0-2's orbit as determined by the fit of 50-2 alone (dash-dotted line) and the 
simultaneous fit of 50-38 and 50-2 (solid line). Note that these twzo best-fit models are very similar on the 
scale of the plots shown here. The middle and bottom rows show the residuals of the position and radial 
velocity measurements from the best-fit models. The middle row shows the residuals from the fit of S0-2 
alone and the bottom rows shows the residuals from the simultaneous fit of S0-2 and S0-38. In all plots, the 
detections of 50-2 are shown in different colors and shapes depending on the method in which they were 
detected: AO detections are purple circles and direct detections from the speckle holography images are teal 
squares. The astrometry of 50-2 is plotted as the offset from the position ої Sgr A*-radio, which is defined 
as the origin of our absolute reference frame. Note that the uncertainty in the prediction of the radial velocity 
of S0-2 around its next time of closest approach in 2018 exceeds the limits of the y axis of the lower right 
plot. The uncertainty in the RV at this time ranges between +1400 km 8 ! relative to the best-fit model. 
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Fig. 8.— Probability distributions of the 7 black hole parameters as determined by an orbital fit of S0-2 
alone (blue dashed lines), a fit of 50-38 alone (black dotted lines), and a combined orbital fit of 50-2 and 
50-38 (red solid lines). All three probability distributions in each panel are shifted by the bias determined 
in the jack knife analysis on the construction of the absolute reference frame (see Appendix (C). This shift 
is also applied in all subsequent probability distribution plots. The black, horizontal error bars show the 
uncertainties in this bias shift. The addition of S0-38 to the orbital fit increases the constraints on the black 
hole parameters, most notably the R, and the mass of Sgr A*. 
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Fig. 9.— Joint probability distribution of x, and y, (left), the position of Sgr А" on the plane of the sky, and 
the joint probability of У, and V,, the velocity of Sgr А" on the plane of the sky (right). The joint probability 
distributions as determined by three orbit fits are shown: 50-2 alone (blue, dashed line), S0-38 alone (black, 
dotted line), and S0-2 and 50-38 together (red, solid line). All the probability distributions are shifted by 
the bias determined in the jack knife analysis (see Appendix ІС). The black, solid lines show the 1, 2, and 
3-sigma uncertainties in this bias shift. 


4.2.2. Values Derived from the Simultaneous Orbital Fit of 50-38 and 50-2 


We fit the orbits of S0-38 and S0-2 simultaneously to constrain the gravitational potential. Using the 
information from both of these stars gives the final values of Mp; and Ro. Figure [6| and Figure [7] (top and 
bottom rows) show the resulting best-fit model orbits for 50-38 and 50-2 respectively in AR.A., ADecl., 
and radial velocity versus time. The astrometric data are shown in an absolute reference frame that defines 
Sgr A*-radio to be at rest at the origin. These figures also give the residuals of these data from the model. 
The residuals show that orbital models fit the data well within the uncertainties, with a reduced chi-squared 
of 1.3. Table [4] shows the best fit and errors on the black hole and orbital parameters of S0-2 and S0-38 as 
determined by the simultaneous fit of these stars. 


The best-fit black hole parameters presented in Table [4] from the fits of S0-2 alone, S0-38 alone, and 
S0-2 and S0-38 simultaneously are taken as the weighted average of the MULTINEST samples for each 
parameter. The fitting error is then taken as the weighted standard deviation of the MULTINEST samples. 
The best-fit values are then shifted by the bias determined by a jack knife analysis on the 7 SiO masers used 
to construct the absolute reference frame (see Appendix |C] for details). The error in this bias shift is also 
presented as a second error term in Table |4| All probability distributions of the 7 black hole parameters 
presented in this work include the bias shift. Note that the jack knife uncertainties due to the absolute 
reference construction are small compared to the fitting errors for Mp, and Ко, so systematics in the reference 
frame do not greatly affect the values of these two key black hole parameters (see Figure 8). Additionally, 
the black hole position and velocity on the sky (xo, yo, Vx, and Vy) as determined by the fit of S0-2 and 50-38 
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are offset from 0.0, but they become consistent with 0.0 when the reference frame errors are included (see 
Figure 9). 


The addition of S0-38 leads to a better constraint of the black hole parameters compared to using the 
information from S0-2 alone. The best-fit black hole parameters and the errors on these values for the 
simultaneous fit of 50-2 and 50-38 as well as the fit of 50-2 alone are shown in Table[4] The errors on the 
mass of Sgr A* and the distance to the Galactic Center decrease by a factor of —2 and 2.5 respectively. 
The best fit values of Mpn and R, as determined from the combined 50-2 and 50-38 orbital fit are: Mp, = 
4.02+0.16+0.04 x 109М., and 7.86 +0.14 + 0.04 kpc. Figure[10|shows the joint probability distributions 
of Мьһ and К, as determined by S0-2 alone, 50-38 alone, and by 50-2 and 50-38 together. 
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Fig. 10.— Left: 2D joint probability distribution with 1-, 2-, and 3-sigma contours of Mp, and R, as derived 
by the orbital fit of S0-2 alone (blue dashed lines), 50-38 alone (black dotted line), and a simultaneous fit 
of S0-2 and 50-38, with the new speckle holography detections of 50-38 included (red solid lines). The 
orientation of S0-38's orbit on the plane of the sky leads to a different correlation of Mp, and Ro, resulting 
in small errors on these parameters when S0-2 and S0-38 are fit simultaneously. The precision of our 
measurements of Мьһ and R, has increased by a factor of ~2 and 2.5 respectively compared to S0-2 alone. 
Note that the у, parameter (not visualized here) is also highly correlated with Mp, and Rọ (see Section [.1). 
Right: Final joint probability distribution of Mz, and А, as derived by the simultaneous fit of 50-2 and 
50-38. 


To confirm that the estimated uncertainties on Mpp and R, from the simultaneous fit of 50-2 and 50-38 
are reasonable given the measurement uncertainties, we compare our estimated uncertainties to the distri- 
butions of the uncertainties on Mpp and Ro derived from the 100 mock data sets. These mock data sets are 
simulated assuming the measurement uncertainties are only statistical. The estimated uncertainties оп Mpn 
and R, from the real data are well within the distributions of uncertainties derived from the simulated data 
sets, with values smaller on average than the widths resulting from the mock data by only 1- and 0.7-sigma 
respectively. In addition, the R, uncertainty derived from the real data falls in the most likely bin of the sim- 


аі 


ulated uncertainties. We therefore conclude that the estimated uncertainties on Mpp апа R, аге reasonable 
given the measurement uncertainties. 
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Fig. 11.— Histogram of standard deviation of posterior samples of Mpp and А, generated by the MULTI- 
NEST fit of 100 fake sets of 50-2 and 50-38 data. The standard deviation quantifies the width of the Mp, 
and К, probability distributions. The 100 data sets were created by assuming only statistical measurement 
uncertainties. The vertical dashed lines show standard deviations of Мьһ and К, samples generated by the 
fit of the real 50-2 and 50-38 data, which are smaller on average than the simulated standard deviations by 
1- and 0.7-sigma respectively. 


The mass-to-distance ratio of Sgr A*, Ми) / Ко» sets the apparent size of the black hole's gravitational 
radius at the distance of К, to Sgr A*: rg = GMpn Fi c^R,. The probability distribution of this value as deter- 
mined by Фе orbital fit of 50-2 and 50-38 is shown in Figure[12] We find that the best fit angular size of the 
gravitational radius of Sgr А" is 5.20 +0.12 +0.06, where the second error term is the estimated uncertainty 


in the reference frame, jack knife bias shifts of Мь and R, propagated to the gravitational radius. 


The additional information found in S0-38's orbital motion also gives increased constraints on extended 
dark mass within the apoapse distances of its and S0-2's orbits. In addition to the central supermassive 
black hole, a central cluster of stellar mass black holes and/or other compact objects has been predicted and 
theoretically explored (e.g. | Morrs|1993} Miralda-Fscudé & Gould 000) Freitag et al 2006 MerriiDTO) 
The amount of extended dark mass within the orbits of 50-2 and S0-38 influences their orbital motion. 
Figure [I3] shows the probability distribution of the extended mass contained within this radius for a fit of 
50-2 alone and a simultaneous fit of 50-2 and 50-38. The resulting 68.3% confidence (10) upper limits on 
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Fig. 12.— Probability distribution of the angular size of the gravitational radius of Sgr АЖ: GMpp,/c7Ro. The 
distribution is calculated from the Мьһ and Ro distributions that are shifted by the negligible bias determined 
in the jack knife analysis (see Appendix С). 


the Ma; within 0.011 pc are 0.06 x 10°M. for ће 50-2 only fit and 0.05 x 10°M for the simultaneous S0-2 
апа 50-38 fit. The 99.7% confidence (3c) upper limits for these fits are 0.16 and 0.15 x 106M. respectively. 


5. Discussion 
5.1. Improvements with the Additional Information from S0-38 


We have demonstrated the power of our new methodology of searching for speckle holography detec- 
tions of 50-38 using the information in the AO detections by applying it to 50-38. The combination of S0-2 
and 50-38 leads to a significant improvement in our knowledge of К, and the mass of Sgr А" compared to 
using the information from S0-2 alone. The reasons why the addition of S0-38 improves the constraints on 
Mp, and R, so much are described here. 


The results of fitting S0-38 alone show that S0-38 primarily gives information about the position and 
velocity of Sgr А" on the plane of the sky (Хо, yo, Ух, and №). Figure shows the joint probability 
distributions of these four parameters with Mp, and А, from the fit of 50-2 alone and in the simultaneous 
fit of 50-2 and 50-38. In the fit of S0-2 alone, Mz, and R, are both correlated with yo, demonstrating that a 
better constraint on yọ would lead to a better constraint on Му; and Ro. In the simultaneous fit, the errors on 
yo and thereby Мьһ and А, are significantly reduced. So it is specifically through improved knowledge of 
yo that 50-38 decreases the errors on Mp, and Ro when simultaneously fit with SO-2. The fact that Mp, and 
R, are not correlated with V, in Figure[14] additionally shows that the difference in values of this parameter 
from the separate fits of 50-2 and 50-38 does not affect the Mpp and Ro results from the simultaneous fit. 


The reason that the addition of 50-38 gives such an improved constraint on y, is the orientation of this 
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Fig. 13.— The probability distribution of the extended dark mass within 0.011 pc of the supermassive black 
hole, as determined by the fit of SO-2 alone (blue dotted line) and by the simultaneous fit of S0-2 and 50-38 
(red solid line). The 99.7% confidence upper limit of the extended mass component decreases by ~10% 
when the information contained in S0-38's orbital motion is added. 


star's orbit on the plane of the sky. Due to S0-2's orientation on the plane of the sky and the fact that we must 
omit the detections of 50-2 around its periapse position because of confusion with the МІК counterpart of Sgr 
A*, S0-2 gives a better constraint on x, than yo. Figure[5]shows the current set of astrometric measurements 
for S0-2, including those points left out of the fit due to confusion with Sgr A* and other known sources 
(indicated by open circles), leaving very few unbiased points in the lower third of the ellipse of SO-2’s orbit. 
50-38 is fortunately on an orbit that is nearly perpendicular on the plane of the sky to that of 50-2, also 
shown in Figure[5] Therefore, this star can further constrain the у, position of Sgr А" and in turn the values 
of Mp, and Ro. The best-fit solutions of Ми, and R, for each case agree well within 1 sigma. With the 
addition of the information provided by 50-38, the errors on these parameters decrease by a factor of ~2 
and 2.5 respectively. 


The joint probability distribution of Mp, and Ro (Figure [10) also shows that S0-38 contains different 
information than 50-2 about these parameters. Mp, and К, are not individually constrained by the orbital 
fit of 50-38 alone, but their correlation is constrained. The correlation as determined by 50-38 alone has a 
different slope than the correlation determined by S0-2 alone. This different slope leads to the additional 
constraints on Mpp and Rọ when information from 50-38 is added. In the future, the addition of more stars 
with high orbital phase coverage and orbits perpendicular to S0-2 will lead to even smaller errors on Mpp 
and Ro. 


L4. 


5.2. Scientific Implications of New Constraints оп М, and R, 


The values of Mpp and А, presented here agree within uncertainties with previous measurements made 
using stellar orbits. The most recent previous measurements from Keck and VLT data respectively are 


4.1+0.4 (Meyer et al.|2012) and 4.30 + 0.36 (Genzel et al.|2010) million solar masses for the mass of 
the black hole and 7.7 + 0.4 (Meyer et al.|2012) and 8.2 + 0.34 (Gillessen et al.|2013) kpc for Ro. The 


uncertainties derived from the orbital fit of S0-2 alone are slightly smaller than the uncertainties of these 


previous measurements due to the increased orbital phase coverage of this star, but the primary improvement 
in this work is the added information from the orbit ої 50-38. 


The method of using multiple stars’ orbits to determine Мьһ and R, was also used in|Gillessen et al.] 
(2009b). In that work, the orbits of 5 stars in addition to S0-2 were simultaneously fit to determine the 
gravitational potential parameters. These 5 stars had orbital phase coverages ranging from 12 - 35% and 
three of the stars had multiple radial velocity measurements. S0-38 was not included in this set of stars, 
since it had only been observed from 2005 through 2008 at the time. With our observations now covering 
over 80% of S0-38's orbit, just including this one star in addition to 50-2 has significantly decreased the 
errors on Мьһ and Ro. From the arguments presented in the previous section, this seems to mainly be due 
to S0-38's orientation on the sky, as well as the more than 8096 orbital phase coverage. Of the 5 other 


stars used in Gillessen et al.|(2009b), two of them have orientations ~ 45 degrees away from S0-2's vertical 


orientation, but none has the perpendicular orientation of 50-38. 


The R, value presented in this work also agrees with other, recent direct measurements of the distance 
to the center of our galaxy within ~ 2 sigma. VLBI measurements of the trigonometric parallax of H20 


masers in the star forming region Sgr B2 give a value of 7.9-- 0.8 kpc (Reid et al.]2009), while dynamical 
modelling of the nuclear star cluster gives a statistical parallax distance to the Galactic Center of 8.27 +0.13 


kpc (Chatzopoulos et al.|2015). Another recent indirect measurement, using VLBI parallaxes and proper 


motions of over 100 masers to model the motion and structure of the Milky Way, gives a comparable value 
and uncertainty in А, to that of the statistical parallax method: 8.34 + 0.16 Крс (Reid et al.|2014). The 
direct statistical parallax method and indirect modelling method have similar uncertainties in R,, but a 


value of ~ 2 sigma higher than presented in this work. It remains to be seen whether continued orbital 
monitoring and further improvement of К, constraints from the orbits ої S0-2 and S0-38 will maintain this 
mild disagreement. 


The direct measurement of К, presented here has implications for constraints on the structure and 
kinematics of our galaxy. measured the apparent proper motion of Sgr A*, which 
is due only to the galactic orbit of the sun if the supermassive black hole is assumed to be at the center of 
the galaxy. This measurement gave a ratio of the circular rotation speed at the radius of the Sun (O,) and Ro 
of 29.45 +0.15 km s! кре !. Combining this ratio with the R, measured in this work gives Ө, = 231 +4.3 


km s^!. This value of ©, agrees well with the independently measured value found by 2014| of 
240 +8 km 8. 


The new mass value presented here does not significantly change the position of the Milky Way in 
the observed correlations between mass of the central black hole and host galaxy properties, such as stellar 
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velocity dispersion and mass of the bulge (see|Kormendy & Ho|2013|for a review). Sgr А", along with other 


central black holes in galaxies with pseudobulges, has a lower mass than expected from the tight correlation 
seen in ellipticals and galaxies with bulges. 


The new upper limit on the extended dark mass within the orbits of S0-2 and S0-38 is also improved 
compared to previous work. To compare the limits presented in Section 4] with previous measurements from 
and (2009b), we tranform our upper limits to find the extended mass 
within the apoapse distance of 50-2: 3.1 x 10!! km = 0.01 pc. The 1-sigma (3-sigma) upper limits within 
this radius are 0.05 (0.14) x 10°Mo for the S0-2 only fit and 0.04 (0.13) x 109М for the simultaneous 50-2 
and 50-38 fit. and {Gillessen et al.|(2009b) found 1-sigma upper limits of 0.12 and 0.17 
x 106M. respectively; therefore our new upper limits are a factor of 3-4 lower than previous measurements. 
50-38 does improve the limit compared to S0-2 alone, but the main reason for the lower limit compared to 
earlier work is due to increased time coverage of the orbit of S0-2. Our new limit is still over an order of 
magnitude more than the ~500 - 1000 Mo of stellar remnants predicted to be within 0.01 pc (e.g., 
[et аї |2006), and other models predict even less mass within 0.01 pc (Merritt|2010). As observations of 
the orbits of S0-2 and 50-38 continue, the limits on the extended dark mass within 0.01 pc will continue to 
decrease. 


The improved constraints on the gravitational potential that come with the addition of S0-38 also im- 
pacts future tests of general relativity in the Galactic Center. General relativistic deviations from pure 
Keplerian motion are expected to be detectable when S0-2 goes though its time of closest approach to the 
black hole in 2018. The deviations will be observable as the shift of the measured velocity of S0-2 due to 
the gravitational redshift. Measuring the deviations from SO-2’s Keplerian orbit require as much knowledge 
of SO-2’s Keplerian orbit and the gravitational potential as possible, so the additional constraints from 50- 
38 are important to this future probe of general relativity. Another observable deviation from a Keplerian 
orbit is the precession of the point of periapse. This general relativistic precession is confounded by the 
Newtonian precession due to extended dark mass within stellar orbits. Measuring the general relativistic 
precession therefore requires the measurement of precession in at least two stars. Additionally, the mea- 
surement of the ratio of Мьһ and А, from stellar orbits is required to compare the size of the black hole 
shadow as measured by the upcoming Event Horizon Telescope to theoretical predictions 
[2015). The theoretical half-opening angle of the shadow of Sgr A* as observed from Earth is predicted to be 
(5+0.2)GMpn /Roc’, regardless of the spin of the black hole. With the addition of 50-38 and an increased 
time baseline of observations, the gravitational radius is now known within ~ 3%, which is less than +4% 


range in the theoretically predicted sizes of the black hole shadow assuming no knowledge of the spin of 
Sgr A*. In the future, our knowledge of the gravitational potential in the Galactic Center will increase with 
more observations of S0-2 and S0-38 as well as with the addition of other short-period stars, thereby also 
increasing the possibility of measuring the effects of general relativity in this extreme environment. 
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Fig. 14.— Joint probability distributions of Му, (top) and А, (bottom) and хо, yo, Vx, and V, as determined 
by the orbital fit of S0-2 alone (blue dotted line) and by the simultaneous orbital fit of 50-2 and 50-38, with 
all new speckle holography detections included (red solid line). These probability distributions are shifted 
by the bias determined in the jack knife analysis (see Appendix [C). The correlation between these sets of 
parameters in the case of 50-2 alone shows that our knowledge of Мьһ and Ro is limited by how well S0-2 
constrains yọ. With the addition of the measurements of S0-38's position and КУ, yo, and therefore also Mp, 
and Ro, are much better determined. The errors in Mpp and R, in the combined 50-2 and 50-38 fit are ~2 
and ~2.5 times smaller than in the SO-2-only fit respectively. 
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A. Astrometry ої Secondary Standards 


The astrometric absolute reference frame is updated using the methods described in| Yelda et al.|(2014) 
with the addition of new observations of the masers, which are summarized in Table The results of 


this analysis, carried out by Sylvana Yelda, are the updated astrometry measurements for the IR secondary 
astrometric standards, shown in Table|[6] 


Table 5. Summary of New Maser Mosaic Observations 
Date Date Frames Frames FWHM Strehl Nstars Кіті с PE 
(UT) (Decimal) Obtained Used (mas) (mag) mas 
2011 July 20 2011.549 65 64 62 0.23 2103 15.8 1.87 
2012 May 16 2012.373 54 54 58 0.24 2040 15.7 1.82 
2013 July 1 2013.501 175 173 59 0.24 2685 16.4 1.63 


“Kim is the magnitude at which the cumulative distribution function of the observed K magnitudes reaches 
90% of the total sample size. 


bPositional error taken as error on the mean from the three sub-images in each epoch and is derived from stars 
with K « 15. These positional errors also include distortion error. 
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Table 6. Galactic Center Secondary IR Astrometric Standards 


Name K To.IR Radius A К.А. ora" ADec. ODec? үрд? Урес 
(mag) (year) (arcsec) (arcsec) (mas) (arcsec) (mas) (mas yr!) (mas yr!) 
S0-6 14.2 2010.11 0.36 0.0200 1.1 -0.3558 Ld -5.2 + 0.1 3.5 + 0.2 
50-11 15.4 2010.42 0.50 0.4934 1.1 -0.0607 1.2 -3.8 + 0.2 -2.4 + 0.3 
50-7 15.4 2010.08 0.52 0.5145 1.1 0.1013 1.2 5,8 +0.2 1.0 + 0.3 


Note. — Table[€]is published in its entirety in the electronic version of this paper. 


?Positional errors include centroiding, alignment, and residual distortion (1 mas) errors, but do not include error in 
position of Sgr A* (2.0 mas, 1.4 mas in RA and Dec, respectively). 


Усіосіїу errors do not include error in velocity of Sgr A* (0.13 mas yr“, 0.23 mas уг"! in RA and Dec, respectively). 
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В. 50-2 Data and Orbital Analysis 


Here we report new astrometric and spectroscopic observations of the star 50-2. New AO astrometric 
measurements were obtained in the same way as described in Section [2.2.1] for new observation epochs. 
Table [7] lists the astrometric measurements of S0-2 used in this work. The speckle holography detections 
used in the fit of SO-2 are all from the set of direct detections discussed in Section.1] i.e., S0-2 was detected 
in all three subset images as well as the main image. Speckle holography direct detections in which S0-2 
was confused with another known source are not used in the orbital fit, as in|Ghez et al.](2008). This is the 
first paper in which we fit SO-2’s orbit with the speckle holography astrometric data. These data supercede 
the shift-and-add speckle data reported in previous works. Table [7] lists the astrometric measurements of 
S0-2 used in this work. 


New spectroscopic measurements of S0-2 were obtained with OSIRIS on the Keck 1 and Keck 2 
telescopes using the Kn3 filter that is centered on the Вг-у line at 2.1661 џт. The previously unreported 
measurements as well as the derived LSR-corrected radial velocities of SO-2 are summarized in Table 
Orbital fits of S0-2 performed in this work used these new radial velocities as well as previously reported 


radial velocities from both Keck and VLT (originally published in 2008| and Gillessen et al. 
2009b} also see|Gillessen et al.|2009a|in which both data sets are presented together). 


In order to perform the orbital fits of stars in our sample (as detailed in Section [3.3.1]. S0-2's data must 
first be fit alone in order to determine the probability distribution of the 7 BH parameters (in addition to 
S0-2's 6 own Keplerian orbital parameters). This 13-dimensional orbital fit was done using the Bayesian 
multimodal nested sampling algorithm called MULTINEST (see and 
2009). We also use the results of this orbital fit of S0-2 alone to compare with the results of fitting 50-38 
alone and 50-2 and 50-38 simultaneously (Sectionß.4). In addition to a purely Keplerian orbital fit, we also 
fit SO-2’s orbit alone with an additional free parameter describing the amount of extended dark mass within 
a characteristic radius. This orbital fit also serves as a comparison to the results from fitting S0-2 and S0-38 
simultaneously. 
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Table 7. 80-2 Astrometric Measurements 
Date AR.A. ADec. AR.A. Error ADec. Error 

(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 
1995.439 -0.0401 0.16216 0.0010 0.00086 
1996.485 -0.0481 0.1553 0.0044 0.0029 
1997.367 -0.0547 0.1388 0.0012 0.0015 
1999.333 -0.0682 0.09808 0.0011 0.00074 
1999.559 -0.0677 0.08971 0.0011 0.00077 
2000.305 -0.0638 0.0747 0.0013 0.0016 
2000.381 -0.06411 0.0643 0.00087 0.0019 
2000.548 -0.0631 0.0563 0.0014 0.0023 
2000.797 -0.0584 0.0519 0.0030 0.0019 
2001.351 -0.0520 0.0261 0.0015 0.0013 
2001.572 -0.0504 0.0152 0.0010 0.0012 
2003.303 0.0386 0.0726 0.0017 0.0016 
2003.554 0.03863 0.08430 0.00089 0.00082 
2003.682 0.0435 0.0949 0.0024 0.0019 
2004.327 0.03595 0.11539 0.00076 0.00062 
2004.564 0.03261 0.12328 0.00079 0.00053 
2004.567 0.0355 0.1261 0.0050 0.0046 
2004.660 0.03137 0.12584 0.00069 0.00060 
2005.312 0.02483 0.14110 0.00081 0.00059 
2005.495 0.0201 0.1473 0.0053 0.0023 
2005.566 0.02084 0.1492 0.00083 0.0011 
2005.580 0.0235 0.1493 0.0038 0.0010 
2006.336 0.01296 0.16191 0.00015 0.00016 
2008.371 -0.01016 | 0.18101 0.00013 0.00014 
2008.562 -0.01232 0.18184 0.00016 0.00016 
2009.340 -0.02114 | 0.18266 0.00011 0.00011 
2009.561 -0.02343 0.18269 0.00013 0.00020 
2009.689 -0.02478 0.18265 0.00018 0.00014 
2010.342 -0.03194 0.18059 0.00012 0.00012 
2010.511 -0.03382 0.17984 0.00013 0.00012 
2010.620 -0.03501 0.17931 0.00016 0.00013 
2011.401 -0.04277 0.17377 0.00022 0.00017 
2011.543 -0.04423 0.17196 0.00014 0.00014 
2011.642 -0.04517 0.17121 0.00025 0.00025 
2012.371 -0.05136 0.16326 0.00013 0.00017 
2012.562 -0.05312 0.16074 0.00014 0.00020 
2013.318 -0.05837 0.14965 0.00021 0.00026 
2013.550 -0.05979 0.14541 0.00014 0.00021 
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Table 8. Summary of New Spectroscopic Observations and Radial Velocities of 50-2 


Date Nframes X tin FWHM? A V. Error Visr° 

(UT) (Decimal) (mas) (kms!) (kms!) (kms) 
2008 May 16 2008.373 11 x 900 sec 75 -417 32 26 
2008 July 25 2008.564 10 x 900 sec 84 -380 43 -7 
2009 May 5 & 6 2009.344 24 x 900 sec 74 -285 32 30 
2010 May 5 & 8 2010.349 16x 900 sec 69 -123 22 30 
2011 July 10 2011.521 12x 900 sec 86 23.5 22 -11 
2012 June 8 - 11 2012441 10 x 900 sec 72 183 49 14 
2012 July 21 & 22 2012.556 11 x 900 sec 102 187 33 -5 
2012 August 12 & 13 2012.616 26 x 900 sec 61 167 25 -13 
2013 May 14 & 16 2013.369 18 x 900 sec 67 328 66 27 
2013 July 25 - 27 2013.566 26 x 900 sec 84 339 23 -18 
2013 August 10 - 13 2013.612 31 x 900 sec 78 349 16 -13 


? Average FWHM of S0-2 in the mosaic made of all frames, measured by fitting a two-dimensional Gaus- 
sian to the source. 


PAII radial velocities are corrected for the local standard of rest velocity. 


*The local standard of rest velocity used to convert the observed velocities to the local standard of rest. 
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C. Jack Knife Approximation to Reference Frame Uncertainties 


We expect systematic errors in the orbital analysis of stars at the Galactic Center to be dominated by 
inaccuracies in our distortion solution and unaccounted errors in the construction of the absolute reference 
frame. While the former is potentially a significant source of uncertainty, currently its exact contribution is 
unclear and thus is the subject of future research. We also expect non-propagated statistical uncertainty in the 
reference frame to be a source of significant error. Such issues may arise from rigid model assumptions or 
improper propagation of error correlations. In this case, this issue is the general of several hundred secondary 
astrometric standards that have been treated as if their uncertainties are uncorrelated. This assumption is 
incorrect as the secondary standards are derived based on only 7 primary astrometric standards. If these 
correlated uncertainties are not incorporated into the full Bayesian orbital analysis, then the underlying 
statistical variance could produce a systematic bias in the resultant posterior distributions of the black hole 
parameters. 


We can gain some insight into the effect and magnitude of this bias through various resampling method- 
ologies. We use a jack knife resampling methodology to infer, to first order, the random statistical bias of an 
unbiased estimator by utilizing subsets of the full set of SiO maser positions that are used in the construction 
of the absolute reference frame (see for details on the refernece frame). Seven subsets 
were made by systematically excluding one maser position at a time. These subsets were then used to con- 
struct seven new reference frames and a full orbital analysis ої 50-2 was performed[] The 50-2 astrometric 
measurements in each reference frame is presented in Table] The results of the orbital fits performed in 
the seven subset reference frames are then used to infer the statistical bias of some estimator x by 


2001) 
Xpias © (n- 1)(ху-1 7 Xn) (СІ) 


with a variance of 
n 


2 n-l = 2 

0^ (Xn) € EU 2. (холу Хи-т) (C2) 
where п is the sample size (in this case 7), x, is the estimator derived using the full sample, Хи-т,; is the 
estimator derived by the exclusion of the i-th maser, and x,-ı is the average of these subsets. Table [IO] lists 
this bias and the square root of its variance (с) on the inferred average value of the 7 black hole parameters 
considered in this work. In the black hole parameter probability distribution figures presented in Sections А] 
and B|we overlay 6896, 9590, and 9996 contours over their respective posteriors. For the joint posteriors we 
infer the correlations through a multidimensional form of Equation [C2} 


n 


-1 
O^ (Xy, yp) = wb (Xn Li — Xn 1) (Yn 14 7 Yn 1) (C3) 
і=1,ј= 


^Note that the orbital analysis of 50-2 was performed оп ће speckle shift-and-add astrometry of 50-2 and not ће new speckle 
holography astrometry. We do not expect this to affect the results of the maser jack knife analysis, since the absolute reference 
frame construction is the same in the jack knife analysis and in the speckle holography analysis presented in this work. 
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As noted in Section 4.2.2] the deviations of the xo, yo, Ух» Vy posteriors appear to be explained by 
previously unaccounted for reference frame uncertainties, while the Ми, and А, posteriors are not affected 
by these systematics. 


We emphases that the bias estimated by this analysis is only relevant assuming that these uncertainties 
are not already incorporated into our Bayesian analysis. The bias define in Equation[C1]is interpreted as, to 
first order, the difference between a statistic inferred from a hypothetical set of infinite masers and the same 
statistic inferred from our set of seven masers. Ideally, an uncertainty of this type would be inherently incor- 
porated into any robust statistical methodology. If, contrary to our assumption, the uncertainties presented 
in this paper do indeed account for reference frame construction uncertainties, then direct incorporation of 
the jack knife results would overestimate the uncertainties and thereby hide additional systematic effects. 
Thus, here we present our Bayesian uncertainties and jack knife result separately in Table А] 


If our assumption holds and if the Bayesian posteriors are parametrized in the sufficient statistic of хи, 
then we may directly derive the resultant ‘total’ probability of a parameter x as: 


Gli, оо) = | diis Phin + Fin) Ро) (C4) 


where Z (x|Xn + pias) is the posterior parameterized in x, and A(Xpias|Xpias, о?2(х,)) is a normal distribution 
whose parametrization is determined by Equations and In the limit of asymptotic normality, the 
net effect of incorporating this uncertainty is to shift the posterior by Xbias and adding the variance o(x,) 
in quadrature with the derived Bayesian uncertainty. All posteriors plotted with the 68%, 95%, and 99% 
contours in Sections Дапа [аге therefore shifted by the bias presented in Table[10] The best-fit values of the 
black hole parameters presented in ТаЫе аге also shifted by this bias and include the standard deviation of 
the bias as an additional uncertainty term. The most significant effect of the bias shifts and the corresponding 
additional uncertainties is in the black hole's position and velocity on the plane of the sky (xo, yo, Ух, and 
V. 
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Table 9.  50-2 Shift-and-Add Astrometric Measurements Used in Jack Knife Analysis 


Date AR.A. ADec. AR.A. Error ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 


All 7 Masers: 


1995.439 -0.0418 0.1637 0.0011 0.0014 
1996.485 -0.0520 0.1557 0.0046 0.0035 
1997.367 -0.0556 0.1384 0.0012 0.0025 
1999.333 -0.0648 0.0923 0.0011 0.0015 
1999.559 -0.0659 0.08949 0.0014 0.00093 
2000.381 -0.06306 | 0.06357 | 0.00097 0.00081 
2000.548 -0.0626 0.0591 0.0012 0.0012 
2000.797 -0.0615 0.0487 0.0024 0.0032 
2001.351 -0.0534 0.0258 0.0015 0.0016 
2001.572 -0.0492 0.0142 0.0012 0.0017 
2003.303 0.0397 0.0703 0.0015 0.0021 
2003.554 0.0410 0.08209 0.0013 0.00067 
2003.682 0.0393 0.08902 0.0018 0.00084 
2004.327 0.03691 0.1140 0.00081 0.0012 
2004.564 0.0333 0.12245 0.0011 0.00093 
2004.567 0.0352 0.1265 0.0050 0.0046 
2004.660 0.03083 0.1255 0.00098 0.0020 
2005.312 0.02412 0.1430 0.00081 0.0011 
2005.495 0.0199 0.1477 0.0053 0.0023 
2005.566 0.0216 0.1493 0.0021 0.0012 
2005.580 0.0233 0.1497 0.0038 0.0010 
2006.336 0.01280 0.16237 0.00013 0.00013 
2008.371 -0.01039 0.18146 | 0.00015 0.00015 
2008.562 -0.01247 0.182277 | 0.00017 0.00017 
2009.340 -0.02141 0.18308 | 0.00016 0.00016 
2009.561 -0.02368 0.18315 0.00016 0.00022 
2009.689 -0.02491 0.18311 0.00018 0.00018 
2010.342 -0.03216 0.18103 0.00015 0.00018 
2010.511 -0.03415 0.18025 0.00021 0.00020 
2010.620 -0.03527 0.17972 0.00019 0.00018 
2011.401 -0.04295 0.17417 0.00024 0.00022 
2011.543 -0.04441 0.17241 0.00017 0.00022 
2011.642 -0.04519 0.17152 0.00026 0.00023 
2012.371 -0.05162 0.16365 0.00016 0.00026 
2012.562 -0.05335 0.16112 0.00018 0.00029 
2013.318 -0.05847 0.15042 0.00024 0.00041 
2013.550 -0.05995 0.14586 0.00018 0.00032 
No IRS 10EE: 

1995.439 -0.0424 0.1639 0.0013 0.0014 
1996.485 -0.0526 0.1558 0.0047 0.0035 
1997.367 -0.0561 0.1385 0.0011 0.0026 


1999.333 -0.0653 0.0925 0.0013 0.0015 
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Table 9— Continued 


Date AR.A. ADec. AR.A. Error — ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 

1999.559 -0.0663 0.08967 0.0014 0.00099 
2000.381 -0.06335 | 0.06360 0.00090 0.00084 
2000.548 -0.0629 0.0594 0.0011 0.0013 
2000.797 -0.0618 0.0489 0.0024 0.0033 
2001.351 -0.0536 0.0261 0.0014 0.0016 
2001.572 -0.0495 0.0145 0.0012 0.0018 
2003.303 0.0395 0.0705 0.0015 0.0022 
2003.554 0.0408 0.08232 0.0013 0.00071 
2003.682 0.0392 0.08920 0.0018 0.00084 
2004.327 0.03679 0.1142 0.00086 0.0011 
2004.564 0.0332 0.12270 0.0011 0.00089 
2004.567 0.0352 0.1268 0.0050 0.0046 
2004.660 0.03072 0.1258 0.00099 0.0020 
2005.312 0.02406 0.1433 0.00083 0.0012 
2005.495 0.0198 0.1480 0.0053 0.0023 
2005.566 0.0215 0.1495 0.0022 0.0012 
2005.580 0.0231 0.1501 0.0038 0.0010 
2006.336 0.01280 0.16265 0.00013 0.00012 
2008.371 -0.01028 0.18174 0.00015 0.00015 
2008.562 -0.01235 0.18255 0.00017 0.00017 
2009.340 -0.02125 0.18337 0.00016 0.00016 
2009.561 -0.02351 0.18342 0.00015 0.00021 
2009.689 -0.02473 0.18340 0.00018 0.00016 
2010.342 -0.03195 0.18130 0.00014 0.00016 
2010.511 -0.03393 0.18053 0.00020 0.00017 
2010.620 -0.03504 0.18004 0.00019 0.00018 
2011.401 -0.04267 0.17447 0.00024 0.00022 
2011.543 -0.04413 0.17276 0.00017 0.00021 
2011.642 -0.04492 0.17183 0.00025 0.00023 
2012.371 -0.05128 0.16391 0.00015 0.00025 
2012.562 -0.05301 0.16143 0.00017 0.00027 
2013.318 -0.05811 0.15079 0.00023 0.00039 
2013.550 -0.05958 0.14618 0.00017 0.00030 
No IRS 123: 

1995.439 -0.0426 0.1612 0.0017 0.0017 
1996.485 -0.0544 0.1538 0.0049 0.0041 
1997.367 -0.0579 0.1367 0.0017 0.0030 
1999.333 -0.0655 0.0904 0.0012 0.0018 
1999.559 -0.0669 0.0880 0.0014 0.0018 
2000.38 1 -0.0651 0.0617 0.0012 0.0019 
2000.548 -0.0640 0.0584 0.0014 0.0015 
2000.797 -0.0639 0.0486 0.0026 0.0041 
2001.351 -0.0552 0.0249 0.0017 0.0020 
2001.572 -0.0505 0.0145 0.0020 0.0021 
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Table 9—Continued 


Date AR.A. ADec. AR.A. Error — ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 

2003.303 0.0385 0.0702 0.0016 0.0023 
2003.554 0.0391 0.0812 0.0014 0.0013 
2003.682 0.0385 0.0880 0.0012 0.0010 
2004.327 0.03550 0.1135 0.00074 0.0012 
2004.564 0.0320 0.1217 0.0012 0.0011 
2004.567 0.0340 0.1261 0.0050 0.0046 
2004.660 0.0308 0.1251 0.0014 0.0023 
2005.312 0.0227 0.1425 0.0012 0.0015 
2005.495 0.0185 0.1473 0.0053 0.0023 
2005.566 0.0198 0.1488 0.0022 0.0020 
2005.580 0.0218 0.1493 0.0039 0.0011 
2006.336 0.01160 0.16187 0.00013 0.00013 
2008.371 -0.01139 0.18136 | 0.00015 0.00016 
2008.562 -0.01355 0.18213 | 0.00016 0.00017 
2009.340 -0.02251 0.18318 | 0.00016 0.00019 
2009.561 -0.02491 0.18320 0.00021 0.00026 
2009.689 -0.02618 0.18310 0.00020 0.00021 
2010.342 -0.03339 0.18122 0.00018 0.00023 
2010.511 -0.03514 0.18023 0.00022 0.00028 
2010.620 -0.03637 0.17981 0.00021 0.00020 
2011.401 -0.04404 0.17436 0.00027 0.00030 
2011.543 -0.04541 0.17245 0.00021 0.00030 
2011.642 -0.04616 0.17169 0.00034 0.00029 
2012.371 -0.05271 0.16394 0.00021 0.00029 
2012.562 -0.05437 0.16143 0.00022 0.00030 
2013.318 -0.05946 0.15118 0.00027 0.00053 
2013.550 -0.06103 0.14638 0.00022 0.00038 
No IRS 15NE: 

1995.439 -0.0407 0.1636 0.0012 0.0014 
1996.485 -0.0509 0.1557 0.0047 0.0035 
1997.367 -0.0546 0.1384 0.0012 0.0025 
1999.333 -0.0641 0.0926 0.0012 0.0015 
1999.559 -0.0652 0.08971 0.0014 0.00091 
2000.381 -0.06237 0.06384 0.00092 0.00083 
2000.548 -0.0619 0.0596 0.0012 0.0012 
2000.797 -0.0608 0.0492 0.0024 0.0032 
2001.351 -0.0527 0.0264 0.0015 0.0016 
2001.572 -0.0487 0.0148 0.0013 0.0017 
2003.303 0.0401 0.0709 0.0015 0.0021 
2003.554 0.0415 0.08275 0.0014 0.00082 
2003.682 0.0397 0.08961 0.0018 0.00088 
2004.327 0.03729 0.1147 0.00081 0.0012 
2004.564 0.0336 0.12320 0.0011 0.00095 
2004.567 0.0356 0.1272 0.0050 0.0046 
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Table 9—Continued 


Date AR.A. ADec. AR.A. Error — ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 
2004.660 0.03120 0.1263 0.00096 0.0020 
2005.312 0.02439 0.1439 0.00086 0.0011 
2005.495 0.0202 0.1486 0.0053 0.0023 
2005.566 0.0219 0.1501 0.0022 0.0012 
2005.580 0.0235 0.1506 0.0038 0.0010 
2006.336 0.01304 0.16336 0.00012 0.00013 
2008.371 -0.01032 0.18266 0.00014 0.00015 
2008.562 -0.01240 0.18346 0.00017 0.00017 
2009.340 -0.02143 0.18436 | 0.00016 0.00016 
2009.561 -0.02369 0.18449 0.00015 0.00021 
2009.689 -0.02495 0.18445 0.00018 0.00016 
2010.342 -0.03224 0.18243 0.00014 0.00016 
2010.511 -0.03424 0.18167 0.00019 0.00016 
2010.620 -0.03542 0.18119 0.00019 0.00018 
2011.401 -0.04312 0.17567 0.00023 0.00022 
2011.543 -0.04465 0.17396 0.00016 0.00020 
2011.642 -0.04542 0.17307 0.00025 0.00024 
2012.371 -0.05186 0.16529 0.00015 0.00023 
2012.562 -0.05365 0.16282 0.00016 0.00025 
2013.318 -0.05887 0.15217 0.00024 0.00037 
2013.550 -0.06035 0.14767 0.00017 0.00029 
No IRS 17: 

1995.439 -0.0414 0.1643 0.0012 0.0014 
1996.485 -0.0516 0.1562 0.0048 0.0035 
1997.367 -0.0553 0.1388 0.0012 0.0025 
1999.333 -0.0647 0.0927 0.0010 0.0014 
1999.559 -0.0657 0.08990 0.0014 0.00095 
2000.381 -0.06292 0.06384 0.00093 0.00081 
2000.548 -0.0624 0.0595 0.0012 0.0012 
2000.797 -0.0614 0.0491 0.0025 0.0034 
2001.351 -0.0532 0.0262 0.0015 0.0016 
2001.572 -0.0491 0.0146 0.0011 0.0017 
2003.303 0.0397 0.0705 0.0014 0.0021 
2003.554 0.0410 0.08229 0.0014 0.00068 
2003.682 0.0394 0.08915 0.0018 0.00082 
2004.327 0.03696 0.1142 0.00079 0.0011 
2004.564 0.0333 0.12264 0.0011 0.00092 
2004.567 0.0353 0.1267 0.0050 0.0046 
2004.660 0.0309 0.1257 0.0010 0.0020 
2005.312 0.02412 0.1432 0.00083 0.0011 
2005.495 0.0199 0.1479 0.0053 0.0023 
2005.566 0.0216 0.1494 0.0021 0.0011 
2005.580 0.0233 0.1500 0.0038 0.0010 
2006.336 0.01275 0.16250 0.00013 0.00013 
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Table 9— Continued 


Date AR.A. ADec. AR.A. Error — ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 
2008.371 -0.01054 0.18153 0.00014 0.00015 
2008.562 -0.01263 0.18232 0.00017 0.00017 
2009.340 -0.02161 0.18309 0.00016 0.00016 
2009.561 -0.02389 0.18315 0.00014 0.00022 
2009.689 -0.02513 0.18311 0.00018 0.00016 
2010.342 -0.03243 0.18099 0.00014 0.00017 
2010.511 -0.03441 0.18022 0.00019 0.00018 
2010.620 -0.03555 0.17970 0.00019 0.00018 
2011.401 -0.04326 0.17409 0.00023 0.00023 
2011.543 -0.04474 0.17235 0.00016 0.00020 
2011.642 -0.04551 0.17144 0.00026 0.00023 
2012.371 -0.05199 0.16352 0.00015 0.00024 
2012.562 -0.05373 0.16100 0.00016 0.00026 
2013.318 -0.05886 0.15030 0.00023 0.00038 
2013.550 -0.06038 0.14573 0.00017 0.00030 
No IRS 28: 

1995.439 -0.0413 0.1646 0.0014 0.0014 
1996.485 -0.0516 0.1565 0.0047 0.0035 
1997.367 -0.0552 0.1391 0.0012 0.0024 
1999.333 -0.0645 0.0930 0.0012 0.0014 
1999.559 -0.0656 0.09016 0.0015 0.00091 
2000.381 -0.06277 0.06402 0.00093 0.00078 
2000.548 -0.0623 0.0598 0.0012 0.0012 
2000.797 -0.0613 0.0493 0.0024 0.0034 
2001.351 -0.0531 0.0265 0.0015 0.0016 
2001.572 -0.0490 0.0148 0.0012 0.0017 
2003.303 0.0398 0.0707 0.0014 0.0021 
2003.554 0.0411 0.08252 0.0013 0.00075 
2003.682 0.0394 0.08940 0.0018 0.00084 
2004.327 0.03703 0.1144 0.00081 0.0011 
2004.564 0.0334 0.12284 0.0011 0.00092 
2004.567 0.0354 0.1269 0.0050 0.0046 
2004.660 0.0309 0.1259 0.0010 0.0020 
2005.312 0.02422 0.1434 0.00079 0.0011 
2005.495 0.0200 0.1481 0.0053 0.0023 
2005.566 0.0217 0.1496 0.0021 0.0013 
2005.580 0.0233 0.1502 0.0038 0.0010 
2006.336 0.01286 0.16268 0.00012 0.00013 
2008.371 -0.01040 | 0.18163 0.00015 0.00015 
2008.562 -0.01249 0.18243 0.00017 0.00017 
2009.340 -0.02147 0.18318 0.00016 0.00016 
2009.561 -0.02374 0.18325 0.00015 0.00021 
2009.689 -0.02496 0.18321 0.00018 0.00016 
2010.342 -0.03224 0.18109 0.00014 0.00017 
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Table 9—Continued 


Date AR.A. ADec. AR.A. Error — ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 
2010.511 -0.03423 0.18030 0.00021 0.00018 
2010.620 -0.03540 0.17974 0.00019 0.00018 
2011.401 -0.04307 0.17419 0.00023 0.00021 
2011.543 -0.04455 0.17244 0.00017 0.00023 
2011.642 -0.04532 0.17151 0.00026 0.00023 
2012.371 -0.05177 0.16357 0.00016 0.00024 
2012.562 -0.05350 0.16107 0.00018 0.00028 
2013.318 -0.05867 0.15029 0.00024 0.00039 
2013.550 -0.06016 0.14575 0.00018 0.00031 
No IRS 7: 

1995.439 -0.0399 0.1631 0.0017 0.0021 
1996.485 -0.0498 0.1561 0.0055 0.0044 
1997.367 -0.0545 0.1386 0.0013 0.0028 
1999.333 -0.0634 0.0920 0.0012 0.0016 
1999.559 -0.0644 0.0892 0.0017 0.0012 
2000.381 -0.0619 0.06333 0.0012 0.00089 
2000.548 -0.0617 0.0587 0.0012 0.0013 
2000.797 -0.0621 0.0483 0.0033 0.0044 
2001.351 -0.0526 0.0254 0.0015 0.0016 
2001.572 -0.0482 0.0140 0.0012 0.0018 
2003.303 0.0414 0.0694 0.0016 0.0021 
2003.554 0.0420 0.0812 0.0015 0.0010 
2003.682 0.0411 0.08802 0.0020 0.00086 
2004.327 0.03773 0.1133 0.00086 0.0013 
2004.564 0.0341 0.1213 0.0012 0.0010 
2004.567 0.0363 0.1255 0.0050 0.0046 
2004.660 0.0322 0.1248 0.0014 0.0020 
2005.312 0.0253 0.1419 0.0010 0.0012 
2005.495 0.0209 0.1465 0.0053 0.0023 
2005.566 0.0230 0.1483 0.0022 0.0013 
2005.580 0.0240 0.1485 0.0038 0.0010 
2006.336 0.01372 0.16106 0.00013 0.00013 
2008.371 -0.00944 0.17984 0.00016 0.00015 
2008.562 -0.01155 0.18067 0.00017 0.00017 
2009.340 -0.02057 0.18139 0.00018 0.00016 
2009.561 -0.02281 0.18135 0.00017 0.00021 
2009.689 -0.02406 0.18129 0.00018 0.00016 
2010.342 -0.03134 0.17919 0.00016 0.00018 
2010.511 -0.03339 0.17839 0.00023 0.00018 
2010.620 -0.03447 0.17777 0.00022 0.00018 
2011.401 -0.04205 0.17215 0.00023 0.00023 
2011.543 -0.04357 0.17041 0.00017 0.00020 
2011.642 -0.04430 0.16950 0.00031 0.00028 
2012.371 -0.05081 0.16147 0.00016 0.00028 
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Table 9— Continued 


Date AR.A. ADec. AR.A. Error — ADec. Error 
(Decimal) (arcsec) (arcsec) (arcsec) (arcsec) 
2012.562 -0.05254 0.15899 0.00016 0.00029 
2013.318 -0.05768 0.14818 0.00026 0.00046 
2013.550 -0.05916 0.14347 0.00019 0.00034 
No IRS 9: 

1995.439 -0.0426 0.1637 0.0013 0.0014 
1996.485 -0.0527 0.1556 0.0047 0.0036 
1997.367 -0.0563 0.1383 0.0013 0.0026 
1999.333 -0.0655 0.0923 0.0013 0.0015 
1999.559 -0.0664 0.08948 0.0015 0.00095 
2000.381 -0.06354 0.06341 0.00095 0.00086 
2000.548 -0.0631 0.0592 0.0012 0.0012 
2000.797 -0.0620 0.0488 0.0024 0.0032 
2001.351 -0.0538 0.0260 0.0015 0.0016 
2001.572 -0.0498 0.0144 0.0013 0.0017 
2003.303 0.0393 0.0703 0.0014 0.0021 
2003.554 0.0406 0.08216 0.0014 0.00071 
2003.682 0.0390 0.08905 0.0018 0.00082 
2004.327 0.03666 0.1141 0.00083 0.0012 
2004.564 0.0330 0.12256 0.0011 0.00090 
2004.567 0.0350 0.1266 0.0050 0.0046 
2004.660 0.03062 0.1256 0.00096 0.0020 
2005.312 0.02396 0.1431 0.00078 0.0011 
2005.495 0.0197 0.1478 0.0053 0.0023 
2005.566 0.0214 0.1494 0.0021 0.0012 
2005.580 0.02302 0.14978 0.0038 0.00099 
2006.336 0.01269 0.16255 0.00013 0.00013 
2008.371 -0.01037 0.18165 0.00015 0.00015 
2008.562 -0.01245 0.18246 0.00017 0.00017 
2009.340 -0.02135 0.18327 0.00016 0.00015 
2009.561 -0.02359 0.18336 0.00014 0.00021 
2009.689 -0.02482 0.18333 0.00018 0.00016 
2010.342 -0.03202 0.18127 0.00014 0.00016 
2010.511 -0.03401 0.18048 0.00020 0.00017 
2010.620 -0.03515 0.17995 0.00020 0.00018 
2011.401 -0.04275 0.17444 0.00023 0.00021 
2011.543 -0.04422 0.17271 0.00016 0.00021 
2011.642 -0.04499 0.17179 0.00025 0.00022 
2012.371 -0.05135 0.16390 0.00015 0.00024 
2012.562 -0.05306 0.16143 0.00017 0.00025 
2013.318 -0.05818 0.15075 0.00023 0.00038 
2013.550 -0.05963 0.14620 0.00018 0.00028 
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Table 10. Jack Knife Bias and Variance 


Black Hole Parameter (units) Bias Standard Deviation of Bias (c) 
Distance (kpc) 0.01 0.04 
Mass (109Мо) 0.09 0.04 
X Position of Sgr A* (mas) -0.50 1.90 
Y Position of Sgr A* (mas) -0.73 1.23 
X Velocity (mas/yr) -0.08 0.13 
Y Velocity (mas/yr) -0.17 0.22 
Z Velocity (km/sec) 4.07 4.28 
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